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ABSTRACT 



The spatial distribution of power In a nuclear 
reactor tends to follow the neutron distribution because the 
energy resulting from nuclear fission Is released largely 
at the point where neutrons Induce fission. Leakage In an 
unreflected reactor causes both distributions to be higher 
near the reactor center. The power at the center Is limited 
by the maximum temperature which can be tolerated by core ma- 
terials with the result that outer regions operate at less 
than ultimate capacity. Ideally all regions should be oper- 
ated at the maximum tolerable temperature, which would require 
uniform spatla .1 distribution of pow«r. Thus an Important 
Improvement In a reactor's power output can be achieved by 
making the power distribution more nearly uniform. 

One simple, economical method of achieving uniform 
power consists of compensating for the decrease In power near 
the edge of the core by means of a reflector. If the return 
of neutrons from the reflector to the core can be adjusted to 
compensate for the leakage factor, the power distribution can 
be markedly Improved, This method, applied to a reflected 
spherical homogeneous reactor, is studied in this paper. 

The distribution of neutrons In a reactor Is 
specified reasonably well by the diffusion type equations to 
be discussed In Chapters I and II, Multi-group methods of 
solution are used because of the importance of the dis- 
tribution of fast neutrons, Assum.ing two groups of neutrons 
was considered Inadequate for the Intermediate type of re- 
reactors which are studied. Multi-group calctilatlons are so 
time-consuming that the use of computing machines is essential. 
Chapters I and II develop techniques for use with M. I, T,'i 
digital computer Whirlwind I, 



ABSTRACT 



(continued) 



To simplify the problem, a basic reactor is 
selected with highly enriched fuel, sodium coolant, beryl- 
liam moderator and reflector. The reactor is assiimed to be 
homogeneous and in the steady state. The effect of reactor 
poisoning is not considered. Spherical geometry is assumed. 

It was found that by suitable adjustment of the 
material proportions in both core and reflector the average 
power could be raised to within l6fo of the maximum power. 
Results are presented and discussed in Chapter III. 
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REr’LEC'.ED HOJ OGS;i:OuS REACTOR 
liTHODuCTIOM 

I;eutron leaAace frorr. an unreflected reactor results 
in a lov/er concentration of neutrons tov/ard the reactor edee. 
By virtue of the fact that the pov;er resulting from nuclear 
fission Induced by neutrons is largely liberated at the 
point v/here the fission occurs, the spatial distribution of 
povfer also is depressed near the reactor edge. This imposes 
a limitation on a reactor's power output. The naximun power 
density at any point is limited by temperature considerations. 
Outer recions, v^ith lov/er pov/er density, operate at tem- 
peratures lov;er than their ultimate capacity. Elementary 
considerations indicate that the average pov/er in an 
unreflected spherical reactor is only one third of the 
maxim.um. pov/er produced at the reactor center. Thus, if a 
uniform power distribution could bo achieved, an important 
increase in reactor output could be realized. 

Spatial variations of the proportions of the m.atc- 
rials in the core, together v/ith use of a reflector, is one 
method of i.;;proving pov/er distribution. Another is the use 
of a poison to depress the neutron concentration near the 
reactor center. These methods have the disadvantages of 
Increased complexltj’’ and cost, a simple, economical alter- 
native, the use of neutrons reflected back into the core 
from the reflector to compensate for the pov/er depression in 
the outer core, is investigated here. 
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The diffusion equations to be discussed in 
Chapters I and II specify reasonably well the distribution 
of neutrons in a reactor. Since neutron energy varies 
widely and since the coefficients of the terms In the dif- 
fusion equation depend on neutron energy, some simplifi- 
cation is necessary. One method, applied to thermal re- 
actors, assum.es that neutrons are slowed to thermal energies 
Immediately after birth. The coefficients of the terms in 
the diffusion equation are t.hcn determined at thermal 
energies, and a simple analytical solution la pocsible. 

However, in this problem, the number of neutrons 
v/ith energies greater than thermal must be increased in order 
to supply the reflector with fast neutron.s for moderation and 
reflection back into the core to compensate for leakage. 
Therefore, the distribution of higher energy neutrons must 
be known, and one-group calculations are Inadequate. 

Considering tv.'o groups of neutrons — fast and 
slow — would give some representation of fast neutrons which 
leak from the core to the reflector, where they may be 
moderated and returned to the core as slov.' neutrons. How- 
ever, better solutions are possible if more than two groups 
are considered. 

Mxxltl-group calculations require solutions of dif- 
fusion equations for each energy group of neutrons con- 
sidered and are very laborious and time-consuming. The use 
of a high speed digital computer serves to .minimize both of 
these problems. In Chapters I and II, diffusion equations for 
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four and six energy groups respectively are broken down Into 
a set of simple calculations v/hlch a oomouter can perform to 
solve the equations. The actual programmed set of In- 
structions Is used v;lth M. I, T.*s digital computer, 

* 

Whirlwind I. The basic method Is applicable to other 
digital computers by a suitable transformation of the 
specific set of instructions. Whirlwind I requires less 
tlian a minute to solve the equations , but a period of from 
two to three minutes Is required for the computer to trans- 
form the Input of decimal numbers and Interpreted Instruc- 
tions to Its own binary system. This conversion time can be 
reduced by rewriting the program In the basic 'Whirlwind code. 
The Whirlwind I computer will perform this function auto- 
matically and will cut a new tape containing the converted 
program. 

In this paper a basic reactor composition and con- 
figuration are chosen, and the proportions of the various ma- 
terials of the core and reflector varied In order to de- 
termine their effect on the power distribution. The funda- 
mental reactor type studied was one considered suitable for 
high output mobile reactor application. 

The requirements of this type of reactor led to 
the use of core materials capable of withstanding high 
temperatures. A liquid metal coolant was chosen in order to 
avoid having to pressurize the core. Of the liquid metals 
available, sodium was chosen because of Its superior heat 
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removal characterlatics and because of Its compatibility with 
other materials at high temperature with respect to corrosion. 

Once sodium was chosen as the coolant i stainless 
steel became the best material for cladding in view of its 
high temoerature strength and resistance to corrosion by 
sodium, 

A dispersal fuel element of uranlxim dioxide in a 
matrix of stainless steel was chosen to meet high temperature 

requirements and to minimize radiation damage. Calculations 

* 

for a fuel element of this type have been made. 

Beryllium was chosen as a moderator and coolant 
because of its high moderating ratio and good thermal proper- 
ties. 

The average operating temperature of the reactor 
o 

was taken as 1000 F, which corresponds to a thermal neutron 
energy of 0,0? ev. 

For mathematical simplification, spherical geometry 
was assumed. The reflector thickness was uniform. Both core 
and reflector were assumed to be homogeneous and in the 
steady state. The effect of reactor poisons was not con- 
sidered, Thus an actual reactor would require additional re- 
activity to. overcome the effects of fission product poisons. 



“Proceedings of the International Conference on the Peaceful 
Uses of Atomic Energy," Volume 9» page 561, C, E. Weber and 
H. H. Hlrsch. 
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CHAPTER I 

Method of Finite Differences 
L. N. Hoover 

This chapter will develop a program for solving 
multigroup neutron diffusion problems In spherical geometry 
on the Whirlwind I (WI) digital computer using a finite 
difference approximation to the age diffusion equation. The 
solution will give the power distribution, the neutron flux 
dlstrlbiitlon In each of several energy groups, and the 
critical radius for a reflected spherical nuclear reactor 
having a homogeneous core and a homogeneous reflector. The 
program will be developed for neutrons In four energy groups. 
By suitable modification of the basic program, the number of 
energy groups and the number of regions considered can be 
Increased. Four energy groups represent the maximum number 
of groups which this program will conveniently handle using 
only the Internal storage of WWI. To Increase the number of 
energy groups considered will require using auxiliary 
storage, which has the disadvantage of Increasing the 
machine time required for a solution.^ The foui>-group 
program requires approximately three minutes of computer 
time for a complete solution. For a parameter study, the use 
of a larger number of energy groups was not considered to-be 
advantageous, since very little Improvement In results v/ould 
be achieved by the use of more groups. 



I. Ap;e Diffusion Equation 



The age diffusion equation provides a good de- 
scription of the neutron behavior of a large class of 
nonthermal reactors. The following development is concerned 
with reactor assemblies that are large compared with the 
neutron mean free path, and in which the neutron energy is 
reduced primarily by elastic collisions. The method accounts 
for both fact and thermial fission, for the fission spectrum, 
and for inelastic scattering. 

The basic equation used is the age-diffusion, which 

gives good results for the case of weak capture and a 

2 

cross section that varies slowly with energy, 

-DV^0(u) + 2;(u)0(u) = - -i- S(r,u) (1.1) 

a an 

The term - represents the excess neutrons entering the 

d u 

lethargy interval over those leaving, assuming that the slov/- 
ing down density, q, can be treated as a continuous function 
of the lethargy, u. This excess number of neutrons is 
balanced in the steady state by the leakage of neutrons out of 
the region on interest, v;hlch is ~DV 0(u), v;here 0(u) is the 
flux per unit lethargy interval, and by the neutron absorption, 
2^0 (u) 0(u), The term S(r,u) represents neutrons entering 
lethargy u at a space point r as a result of production by 
fissions, plus inelastic scattering from a higher energy. 

There are advantages to expressing Equation (1.1) 
in terms of the slowing down density, q, rather than the 
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flux, 0 , The relationship between the flxix per unit lethargy 

Interval and the slowing down density in the case of slowly 

2 

varying capture cross section and weak absorption Is 




For energies below the aou?oe energy , the same number of 
neutrons must slow down past each energyi except for the ef- 
fects of leakage and absorption. Hence » for reactors with 
small absorption In the slowing down range, the slowing down 
density tends to be a less sensitive function of energy than 
the flux, 0 , W'hlch v;lll vary more rapidly than q as ^g(u) 
varies* 



In terms of q. Equation (l.l) becomes 






jgjui 









+ 



q(u) = - + 3(r,u) 



In Equation (l.3)» u Is a continuous variable* 

For multigroup calculations. Equation (1*3) is integrated 
over the lethargy range from to U 2 , corresponding to the 
width of the l*th energy group. ^ g(u) Is considered as 
Independent of position and brought in front of the 
Laplaclan operator, The age-dlf fusion equation for 



lethargy group 1 then becomes 



T r ( P I T 7 ( q^ 



avg 



- - Ql 



OUT 



-h q^IN 4. 



|3(r,]^u)dUi 



( 1 .^) 
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where Uj_ = 



tne Group 

t 



U2 

du. The avcraces are v;ith respect to u in 

1-1 

1 j>r 

‘ is the degradation of neutrons into the 

OUT 



i til group from the group of next higlier energy, and is 

the degradation out of the i 'th group. The approximation of 



replacing the averages of products appearing in Equation( 1.4) 

j S(r,u 

/lU 



by the product of averages yields 

2 

^ avg 

avg 



-(A 



. OUT _ IN 

+ qi - qj. 



i avg 
l) du 



1 avg 



(1.5) 



This procedure is justified if either q or the cross sections 
do not change radically over the group. 

To solve Equation (1.5) when and S are known, 

it is necessary to postulate an additional relationship 

IK OUT 

betv/een (q)j_ avg> Qi > Qi • simplest method would 

be to assume a linear relationship for q across the lethargy 
group and to. take an arithmetic average. Hov<’ever, with 
only four Icthargj’’ groups, the large range of lethargy per 
group tends to make the linear relationship inaccurate. A 
suitable method for use v;ith a small number of energy groups 
is to use the relationship 



avg = a.6) 

v/herc f is a number chosen by trial and error for each 

group and region of the assembly so that f times the spatial 
OUT 

average of q^_ gives a correct value for the spatial 
average of (q)j, avg* solution is not sensitive to the 

choices of f, and it is usually possible, on the basis of 
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ex])erlence, to make adequate choices the first time* From 
the problem solution, the values of (q)j_ are plotted 
versus lethargy as In Figure 12, A smooth curve Is faired 
through these values of (q)j_ to represent the actual 
process of continuous slowing down. From this curve, 
refined values of f can be determined and used to Improve the 
solution. If the Initial estimate was greatly In error* 

II, Average Nuclear Properties of Material In a Lethargy 
Group 

Equation (1*5) expresses the age diffusion equation 
for lethargy group 1 In terms of material properties averaged 
over the I'th group. Within each group, the neutrons are 
assumed to diffuse without energy loss until they have suf- 
fered the average number of collisions which would be 
required to decrease their energy to that of the next lov/er 
group. At this ^olnt It Is supposed that the neutrons are 
suddenly transferred to the latter group. This process is 
supposed to continue while the neutron energy Is degraded 
from the group of birth energy to that of lowest (thermal) 
energy* 

Since the energy spectrum of the neutrons In a re- 
actor covers a continuous range from thermal energies to 
about 10 Mev, the group-diffusion concept Is evidently an 
approximation. This la partially overcome by assigning ap- 
propriate average values for each group to the various 
properties, such as the cross section and the diffusion 
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coerricien't of the core and the reflector. 

a. Dif 1 uslon Coeff iclont . The average diff vision 
coefficient, Dj^, for energy group i is derived from a con^ 
sideration of t’ne fast neutron curi^ent density. The 



3 

definition of D.^ is 



°i 



J^^D(E)^0(r.E)dE 






r ^ V 0(r.E)dE 
y Eo 



( 2 . 1 ) 



It can be seen from this equation that has the same value 
at different spavial points v/itliin the region "^elng con- 
sidered only if 0(r,E) is separable into a function of r 
multiplied by a function of E, i.c., if the variation in the 
neutron flux spectrum from point to point in the medium can 
be neglected. This assvirnptlon is usually made for simplicity. 

Hence, with D(E) replaced by i , the expression for 

becomes yE^ , 



Di = 



^1 



0(E)dE 



Eor the asymptotic case v/here E«<xEq, which is assumed to 
apply for* all energy groups considered here using Be as 
the moderator, and for the case of weak absorption 



0(E) 



— aCli— 

Ej2o(E) 



or 0(u) 



q(^3._ 



In lethargy notation v;e find that Equation (2.2) becomes 
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U« ( q ) ( 1 ) 

- av^ “ 0/^tr^^l avg 

^M^T^riiaTg ‘7^’l avg ^''e 

The last step Involves the approximation of replacing the 
averages of products by the product of averages* This is 
justified If either q or the cross sections do not change 
radically over the group* Hence, this Is a criterion wliich 
must be considered In choosing the energy subdivisions of 
the group* Thus the average diffusion coeffioient for 
group 1 is given by 



Da - ^3 /Zgltr^i avg 
" ( 1 ) 

^I^a\ avg 



( 2 * 3 ) 



b* Total Flux Per Energy Group * The total flux 
in energy group 1, NVj_, is defined as 



NV. ~ f0(u)du = ( q(u)du 

V: 3 ® 



Use of the same approximations as for the diffusion coef- 



ficient gives the result, where 

“•/ 

= ^iTOl avg <‘'>1 avg 



du 



i2.k) 



or 
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NVi 

^“*^1 avg - ^T"( l~) (2.^^a) 

avg 

c • Average Abaorptlon Cross Section Per G-roun . 

The average absorption cross section for energy group 1 la 
defined as 

f -2'g^(u)q(u)du 

/■«•'*■ r-fm 

t 

When the saiae approximations as before are made, the average 
absorption cross section for energy group 1 becomes 



_ ^^a'^‘^3^1 avg 



(2.5) 



Spatial Integration of the Group Diffusion Equations . 
Equation (l«5) describes the neutron behavior in 
energy group 1. Substitution In (1.5) of (2.5) plus the 
following relationships developed previously 



„ OUT - ^^^1 avg 
^1 f ^ 

( 1 ) 

■ avg 



avg 

NV^ 

(q)l avg: -p-23 

/■^a^l avg 



( 1 . 6 ) 

(2.3) 

(2.4a) 



gives 
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+ I^Tai ■^ruXil) ! KVi 



(jZs)l avg^ 



- OUT , 

- ‘^1-1 



ys(r,u 

uc 



)du 



(3«1) 



For spherical geometry, the following substitutions will 

/ X 2 

simplify (3*1) Uy eliminating the Laplacian operator V : 



■Fn = 



l^OUT - 



OUT 



^3n = 



ls(rn,u)du 
'Ui 



(3.2a) 

(3.2b) 

(3.2c) 



A new notation is adopted* Let 1 refer to the l*th 
energy group as previously and the subscript n refer to the 
n’th space point located at a distance from the origin of 
the spherical coordinates. As before, the source term ^ 
represents neutrons entering group i at radius r^ as a 
result of inelastic scattering and fissions in all groups* 
Use ( 3 . 2 ) in ( 3 . 1 ) to get 



2 ip 

dr 



la 



1 + 



' / avg 



1^ 

n - On 



^ OUT ^ * 3 ^ 



( 3 . 3 ) 



The space points are chosen to divide the spatial region 
into a set of sublntervals of width 2^r* The finite dif- 



ference approximation 



Jd2 'f, 



dr 



4 






n 



(Z\r)' 



(3.^) 



and Equation (3.3) gives the recurrence relation 




■ .Tf 




- 10 - 




(3.5) 



where 




bi 



(3.5a) 







avg 



(3.5b) 




J 



(3.5c) 



Equation (3.5) can be verified by direct substitution of 
Equation ( 3 .^) into Equation (3.3). The above substitutions 



for the numerical solution of Equation (3.5) are that 

0 at r - 0 and also at the extrapolated end point out- 
side the assembly. Another condition is that the neutron 
flux and cui’rent density must be made continuous at 
boundaries between regions of different material properties* 
Suitable choices of Z\ r must be made for each region so that 
the boundaries will occur at space points. At a boundary 
point, M, located at radius r - r^j between tv;o different 

these boundary conditions give (primes refer to the 
outer region) . 



4 

follow a method developed by Ehrlich and Hur'vitz, 

a. Boundary Conditions * The boundary conditions 



continuity of flux, 



continuity of current density 



where 



J=-p4^=.^(-F+rf) 



It can be shov/n that (see Appendix C) 



(-P + r g) 



■^rn. 



2 At 






- At 



) - 






~ ^ ^ A T^ 

Therefore , the boundary condition for neutron current 
density, J = J' , la 



(3. 



f Fm i <FM+^r + Fh _^r> - F„ _^J.) 

= ^ [f ‘''m+ |(F'.M+/ir'+F’'M-zir'> ~2§Vr 



(F* F ) 

M +Ar M- Ar 



(3.6) 



b. Method qL Solution * The recurrence relation 
for ^Fn» Equation (3*5) » will be solved using a procedure 
suggested by Ehrlich and Hurwltz#^ This procedure is 
particularly convenient for use v;lth the WI digital computer# 
To solve Equation (3.5) subject to the above 
boundary conditions, new variables and /S’ ^ will be intro- 
duced, These variables are related to the homogeneous 
(In = 0) solution of (3*5) » the defining relation- 

ships 






n +* 1 



M±4- and (9 
An n 




When the. above expression for c< is substituted into (3»5)» 
with - 0 , it is found that 
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Tho recurrence relationship for the/9's is 

/Q '^n— 1 

^ ■■ + I 



n o<r 



n 



n 



(3.8) 



which can be verified by direct substitution of the definition 
of 0 ^ into (3.8). A combination of Equations (3*5)> (3.7), 
and (3.8) gives the follov/lng relationship 

- ^n+1 + 



n 



n 



ocT 



(3.9) 



n+1 



This is proved as follov/s; 

Substitute (3.7) and (3.8) for <X and B into (3.9) 
and combine terras. This gives 



_1 

exT 



n 



n 



' ' n-1 

J 



= ^n^^ - Vl 



- I. 



n 



By Equation (3.5), the rigi^t side becomes 
This is then Equation (3.9), when n is replaced by n+1. 

Equation (3.7) is used to solve for the ^X*s, 
starting with ~ since Aq = 0 and therefore, from the 

definition of ex' ^ =00. The next step is to use (3.8) to 

solve for tixe /3 's, starting with = I^. r’inally (3*9) 
is solved for the F 's from the outside to the center of the 
sph.crical reactor asre.nbly , starting v/lth - 0 at the 



extrapolated end point. 

The equations for and /5 have the same raathe- 
matlcal form but different constants in regions of different 
material properties, at the boundary between regions, o' and 
0 are single valued in order to maintain continuity of flux. 
Thus Equation (3.9) is unchanged at a boundary. At the 
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bouadary point M, at radius Pj. , separating two different 
nedia, the boundary condition for continuity of current 
density can be satisfied by using Equation (3»6). The sub- 
stitution of the recuri’ence relationships for and (5 into 
(3»6) leads to the following relationships across boundary 
point K (primes refer to the outer media); 



where 




( 3 . 10 ) 



(3.11) 



(3.12) 



c. Thermal G-roup . The thermal neutrons are as- 
sumed to be in a group of essentially constant lethargy. It 
is assuuied that neutrons enter the therma.1 group only by 
elastic scattering from the group of next higher energy. 

The thermal group is treated the sa;?.e as other groups except 
for changing the definition of (vriiere the subscript T 

stands for the therr’.al group) . Eor the them’al group 



•5*’ for the case considered here v/lth 13 space points in 
the core so that r = R/15» v/here R is the core radius. 
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NVt s 

4* Energy Q-rouna • Energy Group I is chosen from 
10 Mev to 0.1 Mev (u - 0 to u s 4,606). Glasatone^ plots 
the fission neutron energy spectrum i and It can be seen from 
this plot that very little error Is Introduced by assuming 
that all neutrons released In fission are born Into Group I. 
Groups II and III are chosen to give approximately equal 
lethargy widths down to a thermal energy corresponding to 
1000®F. 



Group 

I 

II 

III 

T 



0 w 4 » 6o6 
4,6o6 - 10.81 
10.81 - 18.80 
18.80 



The energy groups will be as follows; 

Lethargy . u Energy , E Lethargy Width 

10^0,1 Mev 4.6o6 

100..0.2 Kev 6.204 

200-0.070 ev , 7.990 

0.070 ey (looo°F) o 

e. Inelastic Scattering . Provision is made for 
considering Inelastic scattering Into an energy group by the 
definition of the term J S(r:^u)du In the age diffusion 
equation. This term Is defined to be composed of fission 
neutrons born into group 1 plus those Inelastically 
■cattered into group 1 from higher energies. Since a 
neutron must possess a minimum of about 0.1 Mev of energy If 
It is to be Involved in an inelastic scattering process,^ 
Inelastic scattering occurs only in Group I. The group 
widths are large so there will be no inelastic scattering 
across groups. This can be seen from plots of the energy 
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dlstributlon of Inelastlcally scattered neutrons as given In 
the Reactor Handbook, Reference (?)• 

IV. Programming for the Whl rl vl nd ^ Computer * 



of performing slmole mathematical operations (addition, sub- 
traction, multiplication, and dlTlslon) with extreme 
rapidity under the control of a set of Instructions, called 
the program, which la stored In the computer. A detailed 
account of WWI computer operation is given In Chapter II, 
Section II, of this thesis. For additional information, 
the reader Is referred to Reference (l) , which completely 
covers WWI operating techniques. 



chapter develops the relationships which are used In 
programming the WWI solution. The stops of the solution are 
outlined below; 



of fission neutrons born into Group I as a res\ilt of fissions 
in all groups. (This assumed distribution will be Improved 
to the correct value by successive Iterations.) 



Whirlwind I (WWI) is a digital computer capable 



a. Numerical Operations • Section III of this 



(l) Start by assuming the spatial distribution 



(2) From the assumed spatial distribution of fis- 



sion neutrons Into Group I, which Is 




to get 



U 



I 




( 3 ) Calculate the following constants for Group I, 
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UGin^ Equations (3»5) and (5.12); 



^1 



= 2 



ai + 






s - 



JV 

D 



avg 



P = 



1 /Ar ' A 



k = 2 -f- 



2 s { Ar 5 r<- 

V M 

bi(Ar)2 



i: = pk -h 



Ar M k ’ 



+ 1 - 



2 Ar 

3r,. 



T I 

(4) Calculate c<^ starting with cx^2“ b and using 

Equation (3«7) * ex: , = k - — The value of k must be 

n+-i o<-p_ 

corr'oet for tho region considered. At the boundary betv/een 
core and reflector, v/here n = 15, boundary condition ( 3 * 10) 

I -1 

is used, = K 7— * . 

16 bAf 

(5) Calculate starting with 



using equation 



(j,o) , “ oC^ U’ 



V/here 



' n 

h =iArhls, 



D, 



n 



for 1 ^ n - 15 



- 0 

li 



for 16 ^ n ^ 25 
At the boundary, n = 15, Equation (3«H) applies, 

Ar' ^14 



^15 = PI15 + (t + i^)h 5 + 



sAr c< 



15 



(6) Solve for starting v/ith ^P25 “ ^ 



wornin, tov/urd the center of the spherical reactor assembly 



» 



I 



I 

I 

I 



I 
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ualng Equation (3»9)» 

If +1 + /®n 

= * 

c< . , 

n +1 

(7) Calculate the inelastic scattering, ^^7 

from Group I into Group II by use of the relationship 

where JC is the probability that a neutron scattered ln~ 
elastically in Group I will enter Group II. 

(8) Repeat successively steps (2) through (6) for 

Groups II, III, and T, Note that there are two changes to 

the procedure used with Group I: First, the input term, 

II 

S^, is composed only of Inelastlcally scattered neutrons 
entering Group II (since all fission neutrons are born in 
Group I), and is zero for Groups III and T; second, the 
calculation of ^In Equation (3»5o), 



1 



^n 




1-1 



n 



-f 




1-1 OUT 

Includes elastic scattering, , out of the group of 

next higher energy. 

(9) When has been calculated for a.11 energy 

I 

groups, calculate a new source term, Sp, of fission neutrons 
for Group I from the relationship 



^ avg ^h» 

where g^yg Is the average value in group 1. Thus, 

Is proportional to the fission power distribution in the 
reactor. 



I 

I 



I 

I 

i 



I 



-18- 



(10) The nev/ input tern, is normalized and 

used as the basis for a nev; iteration, which repeats the 
above procedure starting v;ith step (3)» 

b. Iteration Procedure , The iteration procedure 
has been divided into tv;o steps: (l) relaxing from the 
assamed input of fission neutrons, to the correct 

spatial distribution; and (2) changing the radius, R, of 
the nev; reactor core to obtain criticality. 

Let the superscript m indicate the number of the 
iteration, v/ith m = 0 indicating initially assumed values. 
Let a circumflex indicate values which have been normalized. 
Step (1) corrects the spatial distribution of ^S^^ 



as follows; 



T 

All iterated input terms '^Sj^ into Group I, starting 



from the assumed arc normalized v/ith respect to an 

* I m-l 

arbitrarily chosen core spatial point, N, whore' S-,j“ is 
normalized to unity. A normalizing factor if is defined such 
that 



^ = 



m-l 
lo m 



lo m 

*37-r 



(4.1) 






and 



ij m 
^n 




rn 

n 



v/here h = the chosen spatial point for normalization to unity 
m = the nunber of the Iteration 
^S^'^ = the calculated incut term to Group I at spatial 

li 



point n 






i 

f 



I 
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Sj^ - 1 = source term of previous Iteration normalized 



normalized source team 



to unity at space point N 

To determine when to stop Iterating, the normalized 



Input term, 3^ , Is compared at each spatial point, n, with 



the preceding Iteration (see figure below). It Is 

^ A_ 







normalization 
point I 



0123 



N . . 15 15 

A m— 1 A m 



specified that at eaoh spatial point the ratio 
must differ from unity by no more than a small value » ^ , 

to permit the program to proceed to the next step. Whenever 
the value ^ Is exceeded at any apace point, the program ’*'111 
reiterate. This la equivalent to requiring that the solution 
for the core power distribution must have converged suffi- 
ciently to permit the solution from one iteration to lie en- 
tirely within an envelope determined by ^ around the so- 
lution of the preceding Iteration. 



continuing to the next step Is based upon the ability of WV/I 
computer to choose one course of action for a positive 
number and a different one for a negative number (conditional 
control transfer). Provision Is made for the criterion to 



The criterion for choosing between Iterating or 
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operate, irrespective of v;hother /^n greater or 

less than unity, by the following scheme; 

(a) When is greater than unity: 

3 

1 - ,XL „ -f € <v 0, not acceptable (reiterate) 

q n 
^n 

S ni-1 

1 - ■ A ■ - - j ■ £ ^0, acceptable (proceed to the 
^ m 

next step) 

3 

(b) When - ■ ■ is less than unity; 

3 ” 

s 

-A -1 + € 4 0, not acceptable (reiterate) 

Q ra 
^n 

q n-1 . 

^n ~1 -{'6 / 0, acceptable (proceed; 

m 

3n 

Step (2) of the iteration procedure Involves changing the 

core radius to get criticality# The procedure is as follows: 

Step (l) above has required that a plot of the 
A m 

Innut term 3j^ must reach a stable shape before preceding to 
step (2)* In step (2) the last calculated value of ^ by 
(Equation 4.1) is used to test for over or under criticality, 
A y less than unity means that the number of fission 
neutrons entering energy G-roup I Increases from Iteration to 
iteration and, therefore, that the multiplication factor, k, 
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l8 greater than unity (and vice versa). For the critical 

_ iUnoJ 

condition, k - 1, the^source term remains constant for suc- 
cessive Iterations, and ^ therefore is unity. 

If Is not v/ithin a small value d of unity, the 
program goes throixgh another iteration using a nev; value of 
the core radius, R, chosen to come closer to criticality. 

The new radius for iteration m Is 

R - y - R 
m - ° m-l ro-1 

Thus, the core radius Is increased If k <Cl, and vice versa. 
The criteria for determining whether the radius Is correct 
or not is: 

(a) When is greater than unity; 

6 - 1^0 not acceptable (modify radius 

and reiterate) 

£ + 1 /* 0 acceptable (proceed to the 

next step) 

(b) When Y Is less than unity; 

y- 1 -A 6 <0 not acceptable 

14-6 >0 acceptable 

When the iteration criteria of both step (1) and 
step (2) are satisfied, the program proceeds to the final 
step of calculating and printing the core power distribution, 
the critical radius, and the neutron flux distribution in 
each energy group, 

V . Flow She e t . 



The program for the solution of the reactor problem 




d 



{ 



I 



f 

( 
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has been described in Section IV. A flow sheet v/111 be 
developed to show schematically the operations discussed 
In Section IV. The following conventions will be used: 

A square box standing on one 
edge denotes a subroutine for 
calculating the quantity within 
the box. 

Circles enclosing numbers are 
used to specify corresponding 
parts of the flow sheets. The 
numbers are called “addresses." 
The circles are called connectors 
and are put at those places In a 
flow diagram to which reference 
will be made at some other point. 





A rectangular box denotes that the 
quantity inside Is to be calcu- 
lated or the operations within 
are to be performed. Frequently 
the formula for calculating the 
quantity is given. 




This symbol indicates a conditional 
transfer to a If A > B and to b 
if A - B. 
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Subscrlpts and superscripts will have the follow- 
ing definitions: 

i denotes an energy group* 
n denotes a space point at radius r^« 
m denotes the number of an iteration, where m = 0i> 
denotes values initially assumed* 

N denotes an arbitrary space point, n - N, chosen 
as a basis for normalization procedures. 

-► the symbol a-^b means that the quantity a replaces 
the quantity b in the storage register in v/hich 
b is stored. 
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PLOW SHEET FOR FOUR G-HOUPS 
Main Pror;ram 










I 



I 



I 

I 

I 



XjC o ^ X* 0i"( i Oi. A'C 3 

i^x-ln Pro ~rar. (Continued) 









STOP 
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boundary 

condition 









I 
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b’LOW brlKET i*'Ox-< FOUR OROUPS 





T“ 

© 






-h 



sR p 



boundary condition 






1 



I 

I 

I 



I 

( 

I 

I 

I 

1 

J 



I 
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iip, ii’u, + 

boundary condition 



l-'Ar 



r o 



/s 



14 



sR 



CX 1 r 




n=lC 



•' This forc of com^o from combination of Equations (1.6), 
(2,4a), (3.2), a,nu (3.5c) \iith Ar = R/l5, '<'hrre R = core 
radius. (There are 15 core space points.) 
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FLOV; SliZET FOR FOUR GROUPS 







Memory Asslprnment (Sec W'v'I program followinp;) 


Find * 


Contents 




pi 


Storage of numbers 


used in" numerical operations 


p2 


Temporary storage 




bl 


R, {i^^f 


)lll, (^r‘) 




(X^in^I’ ^^“^iny 


, 6 , 


11 


Used to indicate number of iterations 


cl 


Group I constants: 


D, C. T. f, 




"^jr' ^1 avg» ^ a » 


2a' ' 


c2 


Group II constants 


; ( Same ) 


c3 


Group III constants; (Same) 



A FLAD is a floating address used to locate a particular 
set of storage registers. 




\ 

I 
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Flad 



c4 

si 

s2 



fl 

fZ 

f3 

f4 

zl 

zZ 

z3 

z4 

xl 

rl 

rZ 

r3 

r4 

r5 

r6 

r7 

r8 

r9 * rlO 



Contents 

Group T constants: (Same) 

I o 

Assumed source 

Calculated source term from fissions in all 

groups jS ” ^n” 

i-I 



n 



II, 



n 



1 -n -25 
l^n 425 



III 



Fn 1^ n^25 



F 14n^25 
n 

Operating registor for storage of registers ci 

when group i is being used 

Temporary storage for ^<=>< 

n 



Temporary storage for 

i-lr 

Temporary storage for I 



n 



II. 



Temporary storage for inelastic scatt., ^^n^in 



Group constants subroutine 

subroutine 

n 

^z3 subroutine 
n 

subroutine 

^Fj^ calculation stored in fl 

^^Fn calculation stored in f2 

calculation stored in f3 
T 

Fj^ calculation stored in f4 
T T 

Calculates ( inelastic scattering and 



store .3 in xl 



I 
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Flad 


Contents 


rl2 


subroutine 


rl3» alO 


Iterates for correct spatial distribution of 
flux 


pl4, all 


Changes radius and reiterates until criticality 
is achieved (both rl3 and rl4 will perforin a 
maximum of 50 iterations before printing out 
the answer so as to avoid a waste of computer 


x2 


time should the solution be divergent) 

I o 

Calculates 3^ from the assumed where 


al and yl 


Control for Group I calculations 


a2 and y2 


Control for Group II calculations 


a3 and y3 


Control for Group III calculations 


a4 and y4 


Control for Group T calculations 


sO 


Calculates 1 where N is spree point 10 

“ q m 


39 


Normalise ^ 3j^® 


9l0 


(See rl3) 


all 


(See rl4) 


sl2 


Puts in r5» r6 , r? » r8 in place of 

Prints core power distribution, flux distri- 
butions, critical radius, ^ , and number of 



iterations 



-52- 



milRLWIIID PROGRAM 

fc TAPE 374-223-008 A HOOVER 
(24,6) 



t=p2 



pi, +1,0 


+2.0 


+3. 


0 


+4.0 


+5.0 


+6.0 


+7.0 


+8.0 


+9. 


0 


+10.0 


+11.0 


+12.0 


+13.0 


+14.0 


+15 


.0 


+16.0 


+0.0 


+50.0 


+50.0 


p2 ,+0 . 0 


+0. 


0 


+0.0 


+0.0 


+0 . 0 


+0.0 


bi,+6o.o 


+0. 


001073 


+0.01444 


+0.052472 




+20.0 


+2.0 


+0. 


0 


+0.00 


+0.03 


+0.000219 


11, +0.0 


+0 , 0 


+0. 


0 


+0.0 


cl, +2. 209 




+1.039 


+0.621 


+0. 


68 


+0.02158 


+0.06884 




+0.000106 




+0 . 


0 


+0.2171 


c2,+0.5756 


+0.522 


+1.528 


+1. 


08 


+0.06535 


+0.13231 


+0.000523 


+0.0 


+0. 16119 


C3, 


+1 . l84 


+0.008951 

+0.544 


+0.537 


+1. 66 


+1.05 

+0.12484 


+0.03813 
c4, +1.190 


+0.12968 


+0.000522 

+1.0 


+1.0 


+0.03734 


+0.1282 


+0.0354 


+0.003914 




+0.0 


sl,+1.33 


+1.31 


+1. 


29 


+1.26 


+1.22 


+1.20 


+1.15 

+0.9S 

s2,+0.0 


+1.10 

+1.14 

DITTO 


+1. 

+1. 

28r 


05 

45 

+0.0 


+1.00 

iliiliiil 


+0.95 


+0.92 


fl,+0.0 


DITTO 


48r 


+0.0 


f2,+0.0 


DITTO 


48r| +0.0 


f3,+o.o 


DITTO 


4 81* 


+0.0 


f4,+0.0 


DITTO 


48ri +0.0 


2l ,+0.0 


DITTO 


40r 


+0.0 


z2,+0.0 


DITTO 


46ri +0.0 


S3, +0.0 


DITTO 


48r 


+0.0 


z4,+0.0 


DITTO 


48rj +0 . 0 


xi,+o.o 


DITTO 


48r 


+0.0 


rl,ita54r 




icazl+16 


idvsl+4 


linrzl+8 


lads 1+1 2 


itszl+18 


lcasl+2 


idvzl 


lt3El+20 


icapl+20 


imrpl+2 


idvpl+4 


imrbl+iO 


ldvzl+20 


idvbl 


itszl+22 


icabl 


idvpl+28 


imrbl 


idvpl+28 


imrsl+l8 


idvzl 


iadpl+2 


itszl+24 


icazl+16 


idv z 1+6 


imrsl+10 


iadsl+l4 


its 


zl+26 


imrbl+10 


imrbl+10 


idv z 1+2 


ladpl+2 


itGSl+28 


icabl+iO 


idvpl+4 


idvbl 


iadpl 


lmrsl+28 


idvpl+2 


ItS' 


t 


lcazl+22 


irarzl+24 


its2t 


iC3pl 


idvzl+20 


iadpl 


imrpl+2 


idvpl+4 


imrbl+10 


Idvbl 


iad2t 


iadt 


itGsl+30 


ispO 


r2,ita25r 


icasl+24 

iadzl+24 


itsz2 

itsz2+2+c 


isc4 


icrl3 

lct5r 


icspl 

icspl 


idvz2+c 

imrpl+28 


idvbl 


ldvzl+20 


idvz2+26 


lmrbl+10 


iadzl+30 


ltGZ2+28 


iscO 


icr9 


icspl 


idvz2+28+c 


iadZl+28 


itsz2+30+c 


lct20r 


IspO 


r3,ita37r 




icabl 

Irarsl 


idvpl+28 


imrbl 


idvpl+28 


idvzl 


itst 


mimm 


itGz3 


isc4 


icrl3 


icat 


imrsl+2+c 


its2t 


icaz3+c 


idvz2+c 


iad2t 


itsz3+2+c 




Ictllr 


icazl+22 


Inirt 


lmrsl+28 


its2t 


icapl+28 


lmrbl+10 


idvbl 


Idv 


Sl+20 


imrz3+26 


idvz2+26 


iad2t 


ltBz3+28 


isc5 


IcrlO 


icaz3+28+c 


ldvz2+28+c 


itss3+30+c 


ict33r 


ispO 


r4,ita75r 




icabl 


ldvpl+28 


imrbl 


idvpl+28 


idvzl 


itot 


icazl+36 


imrzl+4o 


idvzl+32 


its2t 


imrz4 


iadxl 


irru't 


itsz3 


isc4 


icrl3 


ica2t 


imrz4+2+c 




I 
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V,Tiir*lv;lnd Program 



ladxl+ 2 +c 


imrt 


its^t 


icaz 3 +c 


idvz 2 +c 


iad4t 


ltsz 3 + 2 +c 




ictljr 


ica 2 t 


imrz4+28 


iadxl +28 


Imrt 


irarzl +22 


its4t 


icabl +10 


imrbl+iO 


idvzl +2 

imrs4+28 


itst 


lcazl +38 


imrzl+^0 


idvzl +34 


its2t 


iadxl +28 


imrt 


its6t 


icabl+10 


ldvpl+10 


idvbl 


lts8t 


icapl 


ldvpl+2 


iadSt 


imr6t 


its6t 


lcabl+10 


imrpl +28 


idvbl 


idvzl+20 


imrz 3+26 


ldvz2+26 


lad^t 


iadCt 


ltsz 3+28 


lsc5 


icrlO 


ica2t 


ina’z4+30+c 


iadxl+ 30 +c 


Imrt 


lts8t 


icaz3+2S+c 


idvz 2 + 28 +c 


lad8t 


itsz3+30+c 


lct65r 


IspO 


r5.itallr 




icapl +32 


itsfl+48 


iscG 


lcr2^ 


icafl+48+c 


iadz3+46+c 


idvz2+46+c 


Itsf 1+46+c 


lcd2 


ictSr 


IspO 

IlfiUlll 

icr24 


r6,itallr 




leap 1+3 2 


it sf 2+48 


isc6 


icaf2+48+c 


ladz3+46+c 


idvz2+46+c 


Itsf 2+4G4-C 


lcd2 


ict5i'’ 


IspO 


r7.1tallr 


icapl+32 


itsf3+^8 

1 


lsc6 


icr24 


lcaf3+48+c 


iadz3+^6+c 


idvz2+46-fc 


itsf3+46-i-c 


icd2 


lct5r 


ispO 


r8,itallr 




Icapl +32 


ltsf4+i^8 


isco 


icr24 


leaf 4+48+c 


lads3+4C+c 


idv22+46+c 


ltsf4+46+c 


icd2 


ictpr 


ispO 


r9,ltal3r 




isc3 


icrl5 


Icafl+c 


imrbl+12 


itsxl+c 


lct3r 


isc5 


IcrlO 


leaf 1 + 30 +c 


imrbl+l4 


itsxl+ 30 +c 


lct9r 


IspO 


rlO,ita6r 


isc2 


icr25 


icapl +32 


Itsxl+c 


lct3r 


ispO 


rl2.ita2lr 


lcapl +28 


idvbl 


itst 


lsc3 icrl5 


icaz4+c 


imrt 


idvpl+c 


itsz4+c 


ict6r 


lsc5 


icrlO 


icabi+10 


Irarpl+c 

> 


iadbl 


its2t 


lcaz4+30+c 


ldv2t 


itsz^+30+< 


ictl3r 


ispO 


rl3,lsc3 


lcrl5 

icas2+c 


itssl+c 


ict2r 


icail 


iadpl 


itsil 


icail+^ 

rl^,lsc3 


iadpl 


ltsll+4 


1 sup 1+34 


icpal 


ispsll 

iM0A-123.1234t 


icrl5 


icas2+c 


itssl+c 


ict2r 


Icat 


Isupl 


idvpl 


iadpl 


imrbl 


itsbl 


IMOA- 123 . 1234c 


lcall+2 


iadpl 


itsil+2 


lea 11+6 


iadpl 


itsil+6 


i sup 1+3 6 


icpal 


lspsl2 


x2,lsc3 


icrl5 

liiiiiim 


icasl+c 


ldvpl+l8 


Imrpl+c 


Itssl+c 


ict2r 


il 


al,iscO 


icr9 


icacl+c 


itszl+c 


lctal+2 


isprl 


lspr2 


lopr3 


lspr5 


lspr9 


yl,lsc3 


icrl5 


lcabl+l8 


imrf 1+c 


ltss2+c 


ictyl+2 


a2,iscO 


icr9 , 


lcac2+c 


itszl+c 


icta2+2 


Iscl 


lcr4 


icacl+4+c 


itszl+32+c 


lcta2+7 


icacl+l6 


itszl+4o 


lsc2 


lcr25 


leaf i+c 


itsz4+c 


icta2+l4 


isprl 


ispr2 


ispr4 


ispr6 


isprlO 


y2,lsc3 


icrl5 


i cab 1+2 


imrf 2+c 


lads 2+c 


itss2+c 


lcty2+2 


a3,iscO 


icr9 , 


lcac3+c 


itszl+c 


icta3+2 


iscl 


lcr4 


icac2+4+c 


ltszl+32+c 

lcr25 


leaf 2+c 


lcta3+7 

ltsz4+c 


icac2+l6 

icta3+l4 


itszl+4o 

isprl 


lsc2 

lspr2 


lspr4 


IsprT 


y3#isc3 


icrl5 


lcabl+4 


irnrf3+c 
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v,Tiirl\»lnd 


Pro(p?ara 






iads2+c 


itss2+c 


icty3+2 


a4,iscO 


lcr9 


icac4+c 


itazl+c 


icta4+2 


iDCl 


icr4 


icac3f4+c 




itszl+32-fc 


lcta4+7 


icac3+l6 

Icta4+l4 


itazl+4o 


icc2 


icr23 

ispr4 


icaf3+c 


ltsz4+c 


laprl 


lGpr'2 


IsprS 


y4,i3c3 


lcrl5 


1 cab 1+6 


lrnrf4+c 


ladG2+c 


ItaaC+c 


lcty4+2 


s8, icapl 


Id vs 2+10 


itst 


s9,icc3 


icrl3 


icat 


ImrsS+c 


itsG2+c 




icts9+2 


siO,isc3 


lcrl5 


icasl+c 


idvs2+c 


lta2t 


icapl 


iadbl+l6 


iau2t 


icprl3 


icspl 


iadbl+lC 


iad2t 


icprl3 


ictslO+2 


cii,lcabl 


+l6 


iGUt 


iadpi 


IcprlA 


icat 


iaupl 


iadbl+l6 


icprl4 


IFOR 




+c 


+25 


Cl2,isc2 


lcr25 


icafl+c 


ibaz4+c 


ict2r 


lGprl2 


lac 2 


icaz4+c 


itcfl+c 


ict'fi* 


icc2 


icaf2+c 


It sz 4+c 


let Hr 


isprl2 


icc2 


lcr25 


icaz4+c 


itsf2+c 


ictl7i- 


1gc2 


leaf 3+c 


ltS34+C 


lct2ir 


lGprl2 


isc2 


lcaz4+c 


j tsf3+c 


ict26r 


lsc2 


icaf4+c 


ltGz4+C 


ict30r 


isprlS 


isc2 


icazA+c 


It3f4+C 


lct35r- 


isc2 

1MOA-P123 


;Icr23 

.I23’i5f 


leaf i+c 


idvfA 


itof 1+c 




ict40r 


isc2 


loaf 2+c 


idvf4 


lt>r2+c 


inoA 


lct4Cr 

lot52r 


1gc2 

lcaf4 


lcaf3+c 


IdvfA 


ltaf3+c 


IMOA 


ito4t 


iDC2 


lcaf4+c 


iuv^ t 


ltcr4+c 


i^DA 


xct60r 


icaii 


IMOA 


lcail+2 


IMOA 


icabl 


i'lOA 
iWOA 
iSTOr 
i START 


isc3 

lct73r 

T x2 


icrl5 


icarj2+c 


im2*pl+i8 


idvpl+c 



-35- 



Table IV of Appendix B gives the Whirlwind input data 
for runs (1) through (10). The Whirlwind program for run (8) 
has been given on the proceeding pages. The storage assignments 
for the input data are as follows (using run (8) as an 
example) : 

Numerical 

FLAP Value Description of quantity 



bl, 60,0 R, assumed core radius (cm) 

O.OOlOYj) (>^^^) T T 

0.01444 

0.052472 (i'Sj.)^ 

20.0 T, rerlector thickness (cm) 

2.0 ^r' = T/lO (cm) 

0,0 ^^"^in.^l » inelastic scattering (core) 

0.0 * i^^®lQ-stic scattering (reflector) 

0,03 ^ .’iteration criteria for 3% accuracy 

0.000219 



cl, 2.209 

1.039 

0.621 

0.68 

0.02158 

0.06c64 

0.000106 

0.0 

0.2171 



D,' diffusion coefficient, core, Oroup I 
D*, diff. coef., reflector 
f, f factor, core 



f factor, reflector 
s) for core 



^ absorption for reflector 

ij, reciprocal of G-roup I lethargy width 



' for reflector 
y for core absorption 



m-, 



c2. 


(See program) 


Same quantities as in cl. 


Group II 


c3, 


(See Program) 


M It M M tt 

f 


Group III 


c4. 


(See program) 


U It tt It M 

9 


Group T (thermal) 
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< 



I 



CHAPTER II 



Nurnprical integration fiethod 



\i, Kennedy 

I, Oerivation of Equations * 

The followin*^ diffuaion equation is assumed to 



aoply: 




= . 3i 



( 1 ) 



lo^g E 



The subscript i refers to a particular ener/^y group. In 
this case (considering six energy groups) i ranges from 
1 to 6, x’here grouo 6 neutrons are of highest energy, 
group 1 of lovrest (thermal). 

For groun 6 liquation (1) is modified to include 
a fission neutron source. All fission neutron i are 
assumed to be born in group 6. 

aquation (1) is to be made dimensionless, con- 
verted to a standard form (to allox: the computer to use 
one standard calculation form for all groups) and integrated. 
It is first multiplied by ( ref/^ref ) , where Z^ref is a 
reference macroscopic cross-section, here taken as the 
thermal transport cross-section of the reflector. 

Equation (1) becomes; 

'■'This derivation is based on a method developed by a. S, 
Thompson (J. Appl . Phys . 22 , 1223) vrith modifications intro- 
duced to include leakage of the fast group (6) and in order 
to exploit computer advantages. The original method v/as 
developed for hand calculations. 



1 



^ref 2 Xabj Zref ^ 

3 2;ref Itr. ^ ^^ei' ^(2) 

^ <^ lo?;eL 

For conveniencf^ the follovinf^ diraenaionless cross-sections 



are defined; 




ai = 


Z ref / iZtrj_ 


bi = 


2JabSj_/ S^ref 


di = 


2’si/ Tref 


e . = 
1 


Ifi/ £ref 



.vith these definitions hciuatxon ( 2 ) becoraes: 

ai ^2,/ -u V ^ ^ ^i ^rcf /i 



— V ^ 0 . - b . 2 ^ r “ 

Zref ^ ^ ^ 



( 3 ) 



b loggL 



Next, to simplify the energy variable the folloving defi- 
nitions are established: 



^ a ^^ fiss /^ i ^ 

^ ^-' f iss /^ th ^ 



c - (lof^g ^fiss/^^th^ 



-1 



Thus , 



Vdto^E = ^ 



ana the operator (-c 3/c>(^ ) nay be substituted for 
the operator (^/^log EJ in Equation O), giving; 




I 
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3 2 ^ ^ 01 - - 

^ref 






(i^) 



Finally, the operator which is given in spherical co- 

ordinates (with radial variations only) by 



V^- %r^ + 7=-yar 



%- 



becomes, upon considering a dimensionless space parameter 
r/R (where R is the total radius, including core, reflector, 



and extrapolation distance): 






V72 = , 



or 



2 2 - 

^ ^ “ c)('/r)^ 

Thus Equation (4) becomes 

a< 



2 









3 Zref 

2 2 

or, since the operator (R ) In iv-space is equivalent to 
the operator (V ) in r/R-space: 



( 5 ) 



3 2 



ii V ^0^ - bi Z^ef h - ^ t 



ref 









( 6 ) 



Upon multiplying Equation (6) through by 3( 
defining 

A r 3( frefR)^ 

Equation (6) becomes: 

- b.AJfi : 

Note that the parameter A depends on the total radius of 



(7) 



the reactor and is thus not originally known. Instead it is 
assumed (based on one-group calculations) and later corrected 




! 
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ir the asjumfd radiu3 is aub- oi- auper-critical . 

The derivative on the ri3;ht hand side of Kciaation 
(7) be approxii.nted: 

JO that: “ ft = - (u) 

The left hand iidc of equation (h') repreJents los^ fron. the 
i’th Gncr,'i;y n;roap by leakage, absox\otion , and jio’ iny do'-'v'n. 
.he rii^ht haxnd side rroreaentJ the ir.pat to th( i'th p;roup 
conJiJtiny of neutronj Jiowed Iron the next hi^!;h(,r p;roup. 

for i=l ^thcrnal 2 ;roupj Lquation ['c^} is modified 

to: “ ^1^ ^1 “ 

since no neutrons arc lost from the thermal f>:roup by slofiny 
dov/n . 

The input term in Equation (8) is modified for 
energy group 6 (fast) to include fission: 

(The fission of thermal neutrons must be vreighted by an 
energy factor since all other fluxes in («) are per unit 
energy interval and must be multiplit d by this em rgy inter- 
val to become the total flux. The equation is sinolified by 
dividing through by the energj^ interval v/ith the result that 
it apoears as a denominator in the- thermal tern.') 

aquations (8), (8aj , and (*^J are, of the form 

- g^)^i = h 



6 



( 9 ) 






-^^1- 



\'here p combines all loss terms except leaka,<^^c , aiid h repre- 
sents an innut term, for example, considerinf^ ener^^y r^roup 1, 
is given by Ib^T^ j and h by {-/\.^d2^2^ from (6aJ. 

If -chis form of hquations \8) , {8a}, and (9) is 

Ll 2 

laultiplied tnrough by a^(r/3}^ and the operator a^^V 

O 

replaced by its equivalent (a ^ 1 d/d(r/it} (r/ j^' d/d(r/uj 

( r / ^ 

the follo\-ing expression (vdth x'= r/rt) results: 

a^_x^ ^(l/x^}d/dx(aj_x^)d/dx - = aj^x^h 

or, i-'ith = aj^x^ d/dx, 

(Di^ - a^x^g^} = a^x^h (lU) 

Lquation (10) is of the type 
(D^ - ax^+g^) Y = Z 
V'hich may be factored into 

(D + f ) (D - f} y = Z 
provided f satisfies 

Df -f f^ - ax^g^ = 0 
This arises from 

(9 + f}(y - f) Y = Z = (D^ - ax^g^) Y 
D^Y + IDY - D(fY} - f^l s D^Y - axV.% 

D^Y + fDY - (Df)Y - f(DY} - f^Y * D^Y - ax^^r^y 
or, (9f}Y + (f^)Y - (ax^g^)Y = 0 

v'hich gives houation ^12} specifying f. 

In Louation (11} let (D - f)Y = U, so that 
(D f}U = Z 



( 11 ) 



(12) 



( 13 ) 




I 



I 



where 



(D - f)Y = U 



(14) 



Integration of Equation (l4) gives 

Y2 exp{ Jq f dx*/ax'^) exp(-j^^f dx’*/ax'*^)U 



dx"« 



ax 



1 1 r 




(15) 



and from Equation (13) 



U = exp(- Z • -^2 l!I1 + C ) ( 16) 

ax‘2 v° ~x*T2^ irrri? V 

yGc 2 

which combine to give (v/riting m - f dx'/ax' ) 

Y - m Ad/m^) r/^'"ni Z dx'Vax**^ Col dx»“/ax"*2 + C, 

- \J° J (1^, 

The constants and €3 are determined by the conditions 



that 



yi ~ 0 at X - 0 and 
Y = 0 at X ~ 1 



which is equivalent to requiring neutron flux symmetry at 
the reactor center, and requiring the flux to vanish at 
the extrapolated radius. 

By the first condition C2 must be zero, Bj»- the 
second C]_ is given by 



Cl - - (l/m^) m Z dx‘/ax‘^ jdx**/ax»*^ 

30 that Equation (1?) becomes 
Y s (l/m^) m Z dx’/ax'^ jdx* '/ax' 



( 18 ) 



Solution of Equation (18) requires prior calcu- 
lation of f from Equation (12) with the initial condition 
that f(0) = 0 from (ll) and (lO), Equation (12) is written 
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in terns of differences: 

ax^d/dx(f) ♦ (f)^ ~ ax^c;-' = 0 

or, v.'ith subscript 1 indicatinp^ spacf point x and subscript 

2 indicetinc; soace ojinc (x + Ax): 

2 

(f 2 - fx)/Ax = s(o; 2 ^X 2 ^ + glXj^) - 
iolving by quadratic formula gives 

f^ = -^2^2'^/^ ^ ^ ■*■ 2fx/Ax 4 F. 2 ^ 2 ^ * 

+ ^ (19) 

Factors in ^.qaations (1^?) and (10) are sumnarized for 
clarity: 



Energy Group 

1 (th; 

2 

3 

L 

5 

6 



ji ( input ) 




^^(loss) 


(- A ^ 


bj^A 




(- A^od )i? 


b^A 


4 (cd2/A^ )A^ 


(~ A ^ Cd, )0^ 


b^A 


4 (Cd^/A^ )A''5 




b,A 


4 (cd^/Af 


(- A? cd^)^^ 

Ad 

6 

If— — ■— 


bjA 


4 (cd^/A^ )A'^ 


V eicj^i . 

L — 


be A 


4 (cd^/A^ )A? 



i=l 



II. Frogranning of Problem for Computer 

V/hirlwind I (^Av'I) is a digital computer canable of 
performing simple mathematical operations (addition, sub- 



-44- 



traction, multiplication, division) with extreme rapidity 
and under the control of a set of Instructions, called the 
program, which may be stored in the computer* WWI normally 
operates with numbers that are all less than unity, but, by 
accepting a loss in computer speed and storage capacity, may 
be operated with unsealed numbers (thereby decreasing 
enormously the difficulty of setting up the program) in the 
Comprehensive System (CS system) , as is done in this problem* 

In the CS mode there are I36O registers available 
for storage of either instructions or numbers* An in- 
struction occupies one register, but a number occupies two* 

Of the two registers containing a number, the first carries 
a part of the significant figures which represent the 
number expressed in exponential form, while the second car- 
ries the exponent and remaining significant figures* For 
example, the number 123*4 may be expressed as *1234 x 10^, 
with *123 stored in the first register, and the exponent 3 
and the remaining significant *0004 stored in the second* 
Actually, WWI operates with binary numbers, but an undej>- 
standlng of this feature is not essential to the unde2>- 
standing of the problem since the computer accepts decimal 
numbers from the programmer and returns decimal number 
answers* 

The 1360 registers which comprise the internal 
storage capacity of WWI in the C3 mode are called the 
magnetic core memory (MCM). External storage is available 



in the form of magnetic druna and tape but acquisition time 
is large in comparison with the time required for the com- 
puter to obtain a number stored in and is therefore 
not used in this problem. 

Any problem then is limited to a set of numbers 
and instructions v;hich occupy not more than 136u registers. 
In this problem these registers are used (starting at the 
first register available, v.'hich is number 32 J as follows: 

Registers 32-20<^: vrorking storage (for numbers 

ivhich are calculated during 
the problem and must be tem- 
porarily stored) 

Registers 21U-469: storage for 20 values of the 

neutron flux in each of the 
six energy groups 

Registers 470-5^9: storage for 2U space points 

(r/h, = *^^5, .iJ, .15 etc.) 

Registers 510-597: cross-section storage 

Registers 598-639: storage of constants other 

than cross-sections 

Registers 640-end: storage for the pi’ogram 

instructions 

Actually the prograrri for this rroblem uses only some I 250 
registers so that 110 are available for later expansion of 
the problem if desired. 



For the purposes of this problem there are sixteen 
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coranands which WV/I obeys In the CS mode. They each affect 
an operation, designated by three letters, and an address. 

The address may consist simply of the number of a register 
which contains, for example, the number which the computer 
Is to obtain from MCM* However, since a programmer cannot 
know in advance exactly how many registers to set aside. 

It is frequently difficult to specify a register number as 
an address. To eliminate this problem, WI Is designed to 
allov/ use of Indefinite addresses, called flads (for 
floating add resses) and represented by a lettei>-number com- 
bination such as al , a2 , bl , b2, etc. Thus, If a number Is 
to be stored at the end of the program it Is tagged with a 
flad and the flad used In place of the actual register 
number, WWI simply counts the registers as the program is 
read in and remembers the actual address of the flad. Then, 
whenever this flad is referred to later WffL substitutes the 
proper address. 

The three-letter portion of the WWI commands gener- 
ally calls for a certain arithmetic operation. All of these 
occur In the Multiple Register Accumulator (MRA), For ex- 
ample, If a number in storage is to be multiplied by another 
stored nuidber, the first Is taken from MCM to the MRA, The 
second Is then obtained, the required product formed In the 
MRA and then transferred to storage. Various operations, as 
explained below, may either preserve or destroy the originally 
stored quantities. 



The sixteen instructions used in this problem are 



the follo’vlng: 

Instruction Meaning 

its al transfer the con- 

tents of MRA to 
register al 



lex al exchange the coi>- 

tents of MRA and 

register al 

Ira al c.lear MRA and ajid 

to it the contents 

of register al 



Ics al clear MRA and jub- 

traot from it the 
contents of al 

lad al ad d the contents 

of al to MRA 



Isu al subtract the con- 

tents of al from 
MRA 



imr al multiply MRA by al 

and i^ound off 



idv al di2;_ide MRA by 

Contents of al 



Ooeratlon 

the number In the MRA Is 
transferred to register 
al; the former contents of 
al are lost; the contents 
of MRA are unchanged 

the contents of the MRA 
become the new contents of 
register al> and vice versa 

the MRA is cleared (and 
its contents lost) and its 
new content is the number 
stored in al; al remains 
unchanged 

as above except that MRA 
acquires the negative of 
al; al is unchanged 

the contents of al » plus 
the original contents of 
MRA are added and the 
result stored in MRA; al 
is unchanged 

the contents of al are 
subtracted from MRA and 
the result stored in MRA; 
al is unchanged 

the number in MRA is mul- 
tiplied by the number in 
al; the result rounded 
off to the maximum number 
of significant figures 
which can be stored (about 
8 decimally) is retained 
in MRA.; al is unchanged 

MRA is divided by con- 
tents of al , with the 
result stored in MRA; al 
Is unchanged 
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The following two instructions provide transfer of 
control to the instruction stored in the address of the 
instruction. W¥I normally operates sequentially from one 
register to the next. These two Instructions circumvent 
this sequential operation. As an example, the last in- 
struction in an iterative problem might contain a transfer 
of control back to the first in order to perform the 
iteration. 

Meaninp: Oneration 

control transfer 



Instjruction 
isp al 



Icp al 



conditional 
control transfer 



take the next instruction 
from register al and pro- 
ceed from there 

as above, i^ the contents 
of MRA are negative; other- 
V7ise , Ignore, and take the 
next instruction in sequence 



The next five instructions have to do v/lth operations 
performed with a counter. In WVn a counter may be set aside 
to allow repetition of an Instruction, or set of instructions, 
a specified number of times. The specified number is the 
criterion, the number of times the operation has been per- 
formed is the index. WI sets aside both criterion and 
index (with index Initially zero), cycles the index once 
each time the operation is repeated until, upon comparison 
the index is found to equal the criterion (or if the index 
exceeds the criterion, as pointed out below) at which time the 
computer proceeds to the next Instruction in sequence. 



I 



I 

I 
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Instrurtlon Meaning: Operation 

Isc J select counter J select counter J for fol 

lowing operations with 
counter 

(selecting a counter J, by Isc J , does not in Itself set 



criterion and index; this must be done by an additional 

step prior to the first use of counter j.) 

Instruction Meaning Operation 

icr m cycle reset set criterion equal to m, 

index 7,ero 



let al c^yole coun;^ 



id n cycle l^ncrease 

led n C 3 "cle decrease 

(Note that the counting Inst 



increase the index register 
by one and compare with 
criterion; if the index is 
less return to register al; 
otherwise take the next in- 
struction in sequence 

increase the index by n 
(this may set index greater 
than criterion, but is use- 
ful in operating with more 
than one counter) 

decrease index by n 

ction, let al , la the only 



counter command that transfers control. It is therefore 



usually the last instruction in an operation which is to be 
repeated. ) 



During operation with a counter may treat 
numbers (or, in special cases, instructions) which are stored 
in sequence if +c is appended to the address of the first 
number to be treated. For example, to place ten numbers 
which are stored consecutively in ten other consecutive 
registers it suffices to select a counter with criterion ten, 







I 



m 
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Index zero, obtain and tranofer the first naraber, with +c 
appended to the address section of the obtain and transfer 
instructi )ns , and transfer control back to the first instruc 
tion (in, say, register b]_) v;ith an ict bl. bach time the 
coinouter notes the +c in the address it sets c equal to 
trice the index of the counter in use and thus transfers 
the ten desired numbers. The computer sets c to twice the 
index bf cause of the fact that the numbers occupy two 
registers. This hov^ever is not the case if a +c is added 
to the address section of an iso corrunand. Here the address 
must be another instruction, ocouoying one register, so 
the computer conveniently sets c equal to the index. 

The final instruction is used in connection vdtn 
subroutines i-ithin the program. Identical operations which 
are to bt performed at different points in th( program may 
be set uo as a subroutine in order to avoid duplication of 
insti'uction stox'age each time the recuired c.elculation is 
used. Then, at each point the subroutine is needed, an isp 
instruction transfers control to the subroutine. Hov/cver, 
to leave the subroutine and return to th(. proper ooint in 
the main program an isp instruction at the end of the sub- 
routine must have the address of the I'egister immediately 
follovv'ing the register vrhich originally transferred control 
to the subroutine. This is accomplished bv the final instruc 



tion to be considered here 




I 
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Instruction Meaning Operation 

Ita al transfer ajidress the address section of the 

Instruction In al Is changed 
to the register Immediately 
after the register con- 
taining the latest lap In- 
struction which has been 
obeyed 

For example, to obtain the square root of a number at several 
points in the main program, a square root subroutine may be 
set up at the end of the program. Then, each time the 
operation is to be performed the number Is placed In the 
MRA, control transferred to the subroutine, the root calcu- 
lated In the subroutine and placed in the MRA, and control 
transferred back to the main program to the register Im- 
mediately following the original transfer point. 

In the subroutine the first instruction would be 
Ita bl , with bl the last register In the subroutine. In bl 
vrould be stored Isp — , Thus the ita instruction wo_uld 
cause transfer of control back to the main program at the 



proper x>olnt. 

III. Proprranplng of Specific Problem . 



The problem to be solved here, derived in Part I 
of this chapter, may be stated as 




where x Is a dimensionless space parameter r/R total 



- 5 ^- 

a iJ a stored dirnensionles j crosa-acction ; 



Z io the source term, repreaentina; fission neutrons 
(for the highest energy group J or slov;ed neutrons 
(for all lo-'er energy groups) ; and 
f is calculated from a finite difference approximation. 
Fundamentally the coraouter proceeds by the fol- 



io’ -ing steps: 



( 1 ) 

( 2 ) 

( 3 ) 

( 4 ; 



( 5 ) 



( 6 ) 

( 7 ) 



(S) 



Calculate fix); 

hivide f(xj by ax^ and integrate from u to 



X by trapezoidal sumriation; 

fx 

Calculate and store exp 



f dx/ax^ 






Calculate source term Z, initiAlly from 
an assumed flux and critical radius; 

Divide Z by ax^ and multiply by the expo- 
nential computed in step (3), preserving the 



stored exponential for later r(-a^e; 



Integrate 



( L / ax^ ) ( exp ^ f dx/ ax^) 



dx ; 



Divide the above integral by the square of 

the exponential (still preserving the expo- 

2 

nenti'’l) and by ax ; 

rtf verse the order of storage of the quantity 



rteterninf-d in step (7) so that the integral 
from 1 to x may be determined (this step 
is necessary since Vd,I counters incre-ase 
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ratlier than decrease following each repetition). 

Step (8) effectively cancels the minus sign v;hlch 
appears in the input tern to each group, vlth the result 
that the calculated neutron flux is positive, as it should be. 

(9) Beverse the order of storage of the integrals 
determined in step ( 8 ) and multiply by the exj)onential term. 

(10) Store the results of step ( 9 ) and proceed to 
the next energy group , calculations , changing f, a, and Z in 
steps (1) through ( 9 )* 

Following the sixth application of these steps the 
new thermal flux has been computed; WWI proceeds then to the 
next step: 

(11) Compute a new critical radius parameter 



where the primed quantities are the newly computed values. 

Step (11) is based on the fact that once the 
proper solution for the flux is determined 0 ^^^ will not 
change v:ith more Iterations, as long as the radius is equal 
to On the other hand, if the radius assumed was sub- 

or super-critical 0 ^j^ will be less than or more than 0 ^j^ 
(assumed) respectively. 

(12) Compare the tvro t'alues of the radius parameter; 
if they differ by more than an acceptable error recalculate, 
based on the new values of flux and radius; if the reactor 
is critical, as evidenced by a specified small variation in 
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print out specified information flux, radius, pov/er 

distribution) and stop. 

conniJter nay be stooped at the end of a 
problem by an iSTOP instruction. Alternatively, if more 
problems are to be considered they may be listed after 
the first problem and the computer stooped v:ith the special 
instruction sp26. when it reaches this command at the end 
of the first problem the conouter automatically reads in 
the remaininh; problems consecutively. In this vray a 
considerable savins; in time may be realized, since manual 
effort lothei’V\?ise required to start each problem) is mini- 
mized and since many of the constants which arc already in 
storage will not chan?;e from problem to problem, thus 
decreasing; the time required for the computer to convert 
the inout decimal numbers into binary numbers. After the 
final problem an ioTUP command may be read in to stop the 
computer, otherwise the operator will stop the computer 
manuall}’- wh( n, upon completion of the last problem, the 
computer, under the sp26 instruction, calls for more infor- 
mation. 

Ho mentioned in part II of this chapter, the 
numerical values of this problem are stored in the first 
aortion of lych. The first instruction is stored in register 
640. D"velooment of the program for this problem v'ill start 
at this point v-ith numb<^rs specified as they appear. 
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Instruction 

al, iscO 

icrppl 

iscl 

icrpp2 

isc2 

icr20 

isc3 

lcr20 

isc4 

icr6 

ica636 

its598 

isc2 

a 3 > ica430+c 
imr638 
its430+c 
icta3 



f2,isc0 

a2,lca510 

imr254+c 

its34+c 

ica512 

imr298+c 

iad34+c 

lts34+c 

Ica5l4 

imr342+c 

iad34+c 

its34+c 

ica51o 

irar 386 +c 

iad34+c 

its34+c 

ica5l8 

imr430+c 

iad34+c 

imr598 

imr470+c 

iinr470+c 

ItsS^-^-c 

lcta2 



Operation 

Sets counters: 0 for iteration over 
core region, 1 over reflector, 2 and 
3 over entire reactor (20 space points) 
and 4 over six energy groups; ppl and 
op2 are joreset parameters ii^hich are 
set at the opi ning of the program in 
order that the rest of the program 
may remain unchanged 



Stores assumed radius parameter in 598; 
636 contains latest calculated A- 

Scales 0 , to any desired center-line 
value (value initially set in 

register 638 ) 



Select’s counter 0 to iterate over core ; 
double-length registers ^10 through 518 
contain fission parameters; 254 through 
468 contain fluxes; input term Z to 
fast energy ?roup is determined as the 
sum of each multiplied by 

( 3 -) in 598 , x'^ (x in 470 through 508) 
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lca57^ 

its6l8 

lca576 

its620 

ica530 

its622 

ica532 

lts624 



Jtores loss parameters {b-*-§ cd/k^ J and 
a (^ref /?tr) ox' f^roao 6 in v/orkin^ rf,^i 
sters 618-624 (this is done in order 
that other group parurif ters may later 
be 1 laced in these sarie vvorking regi- 
sters and the remainder of the program 
used for all ent rgy gi’cmesj 



Ispml 



a4,lsc4 

icttnl 

ispx6 



ml,l3c0 

a5f ica6l8 

lmr470+c 

imr470+c 

its78+c 

ictaS 

Iscl 

a6, ica620 

imr470+pp3+c 

imr470+pp3+c 

its7o+pp3+c 

icta6 

ica620 

iad6l8 

Imr6l4 

Imr470+pp3 

imr470+pp3 

its78+pp3 

iscO 

a7i ica622 

imr470+c 

imr470+c 

imr598 

itsl22+c 

icta7 

iscl 

a3,ica624 

imr470+pp3+c 

imr470+pp3+c 

imr598 

itsl22+pp3+c 

icta8 



causes a jump around the next 3 instruc- 
tions v^hich provide entry for groujjs 1-5 
a4 iu I'lad to label entry point for groups 
1-5; next register jumped unless new 
has been calculated 

jumps to subroutine in x6 v.dien new 
^th been corntuted 

2 

calculates ax for core region 



as above, for reflector region; pp3 is 
set at the optining oi'" the program to 
tv/ice (because of double -length re-gi- 
stei’S/ the estimated value of Rcore/R^^^^; 
thus, in order to speeix'y a reflector^ 
core combination it is only necessary 
to change ppl,2,3, not prograei 

places av r^g: v xiue of ax'^ at core- 

reflector intei'face; 6I/1. contains O.5, 
lised frequently in averaging 



calculates loss tern b ^ +AjcdA^ ) 

times for core 



as above, for reflector 




I 



Icao^ - 
lau5.?-i- 
lmrol4 

imr^7C'+PP3 
7^+1 p 3 
ifir593 

lt3l2C+pp3 
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as above, at core-rof lector interface 


ica 626 

:lts3P 


p 5 -ac'i ,3 z?ro in rep:* store 39, 3 20 
(j2C- initially set to zero and left 

to act as zero so^arce 


? ca73 

iclv470 

id'/470 


c...•le'alal^'‘J f a.' x-.O-j (..mst be calculated 
separately to avoid division 'o-^ zero) 


iad ' 2 ; 

inir7?> 

iops 2 

its 1.68 

1cs78 

ldv470 

lad 163 

its l 68 


z9 is square root subroutine 


1gc5 

icr^Q 

f lcsi 6 c-i-c 
inirl 68 -fc 
imr470 
idv78-hc 

ladi 66 -i-c 
ladl 6 f--tc 
lmr47G 
laQoi /-i-o 

idv470 

idv470 

iad.l24+c 

ladl 22 +c 

imr 8 o+c 

isps 2 

ltsl 70 +c 

icsoO+c 

idv470 

lad. 7 C+c 

itsl7^+c 

icta3 


calculates reinalning 19 values of f 

/ 


lsc 2 

p . , lcai 68 +c 
idv 7 r>fc 
its i. 6 o+c 
ictgl 


calculates 
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Isc2 

bl, lcalG6+c 

iadl68+c 

Imr6l4 

imr470 

itsl22+c 

Ictbl 


forms areas of f dx/ax^ prior to 
trapezoidal intc=;r' .tion x) 


isc2 

b2,ical22+c 

iadl20+c 

itsl22+c 

ictb2 


sums above areas in traof. soidal inte- 
c;ration {120 pr'eviouslv set to zero) 


lsc2 

bi,ical22+c 

iopsl 

itsl22+c 

ictb3 


/x 2 

calculates exp( / f dx/ax ) ; flad 
zl id'^ntifies Jo 
exponent! il subroutine 


isc2 

b5» ica34+c 

imrl22+c 

it33^+c 

lctb5 


source tern z times exponential 


lsc2 

bo, ica 32 +c 

iad34+c 

imr6l4 

imr470 

itsl68+c 

ictb6 


foi^ms areas of (Z times exponential) 
pr'ior to trapezoidal inte .^ration 


ica 626 
it si 66 


sets 166 to zero 


b7,ical68+c 

ladl66+c 

itsl68+c 

ictb7 


su.ms areas, in trspezoidal intep^ration, 
to compute fx „ 

J Z (exp) dx/stx^ 


iscO 

fl,ical68+c 

ldvl22+c 

idvl22+c 


fou core ’"■P'Piion: 

/x 

/ Z (exn) dx/ax'^ 

Jo 


idv470+c 

ldv470+c 

idv6l8 

it3l68+c 

ictfl 


, ,2 2 
{ exp ) ax 



ical68+pp3 

its7^^ 


-59- 

2 

oets aside (in 74) / 4 (exp) dx/ax 

at corr~rrflectox' Jd interface 


iscl 

g3, ical68+pp3+c 

Idvl22+pp3+c 

idvl22+pp3+c 


for reflector I'egion: 

Z (exp) dx/ax^ 


idv470+pp3+c 

idv470+PP3+c 

idv620 

itsl68+pp3+c 

ictg3 


(exp)^ ax^ 


lca 620 

iad6l8 

imr6l4 

iex74 

Idvl22+pp3 

idvl22H-pp3 

Idv470+pp3 

idv470+pp3 

idv74 , 

Itsl68+pp3 


as abovf , st core-rof lector interface 
(usinfT numerator set aside in 74/ 


isc3 

icr20 

icil9 

kl,isc3 

ical68+c 

icdl 

l3C2 

its34+c 

ictkl 


reverses storage of above (over entire 
n-acto'’) in ordei- to inte.grate from 
1 to X 


lca626 

its74 


sets register 74 to sero 


isc2 

bS, ica 36 +c 

lad34+c 

imr6l4 

imr470 

lts7o+c 

ictb8 


forms areas prior to integration 


lca626 

its 76 


sets register 76 to zero 



- 60 - 



isc2 

b 9 ,ica 78 +c 

iad76+c 

its 78 +c 

ictb9 

lsc3 

icr20 

icll9 

k2, i3c3 

ica 76 +c 

icdl 

isc2 

its34+c 

ictk2 

isc2 

k3, ica3^+c 

its 78 +c 

lctk3 

isc2 

cl, lca70+c 
imrl22+c 
its34+c 
ictcl 

isc4 

isppl+c 

pi, Ispxl 

ispx2 

ispx3 

lspx4 

lspx5 

ispa4 

x6, isc2 
lca626 
c2, iad34+c 
ictc2 
itsll6 

isc2 

ica 626 

c4,iad430+c 

lctc4 

idvll6 

isu 508 

imr 

lad 508 
imr 5 98 
its636 



sums above areas in trapezoidal inte- 
gration 



returns above integral 3 to original 
or.^er of storage 




exp} dx/ax'' 



(exp)"- ax" 



dx 



counter 4 selects energy group by jump 
to xl through x5 lv;hich store through 
02 position constants for next con- 
’^uting 0^ through 0i); the isp aL 
instruction resets counter 4 and pro- 
ceeds to the- 02. storage routine 



to determine radius critic^.ity calcu- 
lates an aporoximation dx 

based on the ne’ ly oonoutrd 0p in the 
registers follo’-ing 34 

as above, based on prrviousl^^ computed 
(or assumed) 0^ 

computes ratio of above integrals 
computes "JT - 1 

multiplies { i( - 1) by any dc3irfd factor 
to control raoidity of conversion, plus 1 
coraoutes new ^ ^ 

^2 = ( 1+\T- y/n) 




I 
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Ica3^ 

Isul. 

Iad3^ 

it3l22 


i;Liif,ar k,ytra? oauted y<^lj.e of 
at re xctor center to nox'r iiize^ 
c:d.cuinted 


isc2 

c6,lca3^+c 

idvl22 

irar 63 B 


nuriniii .e c<*1c j luted ^ to corres- 

pond to orifi^iually chos^rl vauuc (in 
r^p;±3ter oldj at centei’ 


its'fSC'+c 

ictc6 


store normalized in as3i,r;ned 

rericiters 


lcao 36 

Idv6l2 

ispz2 


calculate rt^otal from T'^ 


IdvoOO 

inOA+123.^56c 

icuoSCi 

isu5,?3 

idvG 36 

iad6n4 

Icpal 

lca59B 

1su53u 

id V 5 30 

i.idd04 

icpal 


print out ^’total (sample number given 
to specify number of significant 
figures desired) 

determine criticality by comparing 
the per centage change in A 
vrith an acceptable error stored in 
register 604; if error is too large — 
either plus or minus — return to al 
(the original starting point of the 
program) and recalculate; if error is 
acceetuble, continue to next sequence 
of instructions 


iFOn 

+5 

+0 

+20 


specifies format of solution: blocks 
of 20 numbers arranged 5 pei' line , 
uith a tab betv.'een numbers 


icc2 

c,^ica330+o 
lI!OA+i. 2343070 
ictc7 

ijc2 

c3,ic-.30'-+c 

inOA+1.23'^5f 

ictc8 


prints out 2o values of the calculated 
flux for each of 6 energy groups con- 
sidered 



ipc2 

c lcd342+r' 
iilO .+1.2j45f 
ictc ) 

isc2 

clO,ica29'"’-i-c 

inOA+1.234,5f' 

ictclO 

loc2 

ell,ica254+c 

i^nA+1.2345f‘ 

lsc2 

cl2,ica21C+c 
i'40,.+1.2 34lj: 
ictc 12 
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ifccO 

aiO, lca520 

ImrlJC+c 

ito34+c 

ica522 

inu’3''‘3'0+o 

ij,d34+c 

ItcS^+c 

ica52'+ 

iinr342H-c 

iad34+c 

its34+c 

ica526 

inr23o-t-c 

iad34+c 

its3'4-i-c 

ica528 

lmr254+c 

iad34-fc 

its34+c 

ic trr/. 


calculate the conribution to the 
total pov/er distribution from each 
groixp (fission cross-sections — 
weighted by A^/c for 3 4 5 — 

stored in 520,522,521,526,528) 


ica34 
1 bu 36 
lad 34 


linear extrapolation to determine 
pov;er at center 



itsl22 
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Iscj 

Hi? jica34-fo 

idyl22 

lM0H+.t.23’^af 

ictni? 



dl,icai^r+c 

its2lCM-o 

icfcdl 

iCdSfS 

itc6l3 

lc‘i530 

itj620 

it3C22 

ica536 

lts62-l 

IscC 

d6,ica3-Uc 

lmr470-i-c 

imr 'i-70+c 

irnr 538 

]_rar'-,34 

its3^l-+c 

ictdo 

ica3>+PP3 

itsfc02 

iscl 

e2,ica34+ppjl-^ 

lfarV70+pp3-i-c 

imr47u-i'pp^+c 

iir'.r598 

imr 5 36 

i t y 3 'f")'PP 3"^ *- 

ictc2 

ica53d 

lad55t 

imr6l4 

lmi’C02 

imr't f0-fpp3 
Imr470+pp3 
irar398 

itc3^l-+PP'’ 

lspa4 



print out pov^er distribution 
normalised to unity at center 



insert stopninp; cormnand (e.K»» 
i3T0P or sp26, as discussed 
on page 5k) 

subroutine to be used after 
calculation of stores 0^ 

in assigned registers, places the 
parameters of the next lox;er 
energy group in the working 
registers and calculates the 
source term Z for the next 
lov^er group; reenters program 
at entx'y point set aside on 
page 56 



(THt FOLLOlVlIiG PAGE CONTAINS 
R)UK GIMILAH S'iBdOUTIKES EHIGH 
PEHFOuK THE ,3AI-iL EUJ.CTION AFTER 
GkLGULaTIOK of 0^, 0 , 0^,0^) 
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x?,isc2 

62 , lca3^+i 

lts25''^+c 

lctd2 

101582 

its6lb 

ica534 

its620 

ica538 

its622 

ica 5^0 

itc624 

iscO 

O-T ^lca3^+c 
lor 4 1^0+ c 
inr470+c 
irar503 
imr553 
its34+c 
ic udT 
ica34+pp3 

its 602 
led 

e3» lca3^+PP3+c 

imr470+Dp3+c 

imr470+op3-t-c 

lmr598 

lmr56o 

Its34+pp3+c 

icte3 

lca 558 

lad 56 o 

imr6l4 

imr602 

irnr470+pp3 

imr470+pp3 

lmr593 



x3,lsc2 

d3,lca34+c 

its2j3+c 

lctd3 

lca58o 

IteSlB 

lca338 

1 bs620 

lca542 

it3622 

lca544 

its624 

LscO 

d8j lca34+c 

irar470+c 

lnx'^70+c 

lmr598 

irar562 

ito3^H*c 

lctd8 

ica34+pp3 

its602 

iscl 

e4,ica34+pp3+c 

Imr470+pp3+c 

lmr470+op3+c 

lmr593 

lrar564 

its 34-1- pp3+c 

lcte4 

lca562 

lad 564 

lmrol4 

imr602 

Irar470+pp3 

irnr470+pp3 

imr598 

Ite34+pp3 

lspa4 



x4,isc2 
d4,lca34+c 
its342+c 
lctd4 
lca590 
Its6l8 
ica592 
ltj620 
ica5'^6 
lto622 
lca548 
its624 
iecO 

d9,lca34+c 
imr470+c 
lmr470+c 
imr598 
lmr566 
ite34+c 
ictd9 
Ica34+pp3 

ltG602 
iscl 

e5i Ica34+pp3+c 
imr470+pp3+c 
imr470+pp3+c 
ljnr598 
imr568 
Its34+pp3+c 
ictcp 
lca566 
iad568 
imr6l4 
lrar602 
irar470+pp3 
imr470+pp3 
imr598 
itc344-pp3 
lspa4 
FUl) 

FU 2b) 



x5j lsc2 
d5,io'i 3^+c 
its386+c 

ictdp 

lca594 

Its6l8 

lca596 

its020 

lca550 

itc622 

ica552 

itsu24 

iscO 

cl,ica34+c 

inr470+c 

imr470+c 

1^598 

imr570 

its34+c 

ictel 

Ica34+pp3 

its602 

iscl 

e6,ica34+pp3+c 

Imr470+pp3+c 

imr470+pp3+c 

imr598 

imi-’572 

its34+pp3+c 

icte6 

ica570 

iad572 

imrul4 

imr602 

imr470+pp3 

imr470+pp3 

itnr598 

Its34+pp3 

lspa4 



its34+pp3 
lspa4 

'.if exponential subroutine (LSR 
22, square root subroutine (LSR 

1 start at 640 (terminates read-in and instructs computer 

to take first instruction from register 6U0) 
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Inpu'~- Li :i ArT^an ^ srarn ^ 

(All ?v-.n-numboi'3d r*e;;''is ' -rs include cho followinij. 



odd 



narabened onj 


'■') provide d .uble rc, 5 lcuor storase) 


‘^ogisber ( 3 ) 


Stored Data 


32 - 4 iO 


Initially ceo aside as storage durin 
program calculal.ionG by inserting 




zeroes 


430-43 J 


Insert 20 values of assumed thermal 
flux ejcaled t ■ unity at center 


47 c -5‘^3: 


Insor-G 20 cnosecutive space pjin.s 
(x--=i*/d=. 05 , .1-, . 15 , .... 35 >1.0) 


51- 


Fission pai^au^^Oers : 'v’^rti5c 


512 


'• 4 c 


514 


» 3c 


51a 

513 




520 


Fissl :>ii cross-sec^' .'ns: Ic 



522 

52 i^- 

520 

523 

'■^30 

532 

53^ 

r, o •' 

n !- 

< V.. 

5^0./ 

542 

544 

5 'iG 

5 it 2 



550 

CPO 



554 
55o 
55 ^ 
50 u 
3.^2 
5o4 
5 O A 
500 

5Y0 

572 



'V oj.-iiv dlCui'lL'., 



ui. 



/c 2 c 

'jC 

4 c 
3 c 



i.. )3 3 0 :;r:..nic oon 

( :.acl’idJ 3 1 jCe 
Zv jup 1 0 > di 
tv’oiorp''. i n me; 
si v.vin ' d .;n 
ap, excep:. 
nerinal ;:r 'Uo) 



r.s: 






I* 

h 

II 

li 

II 

It 

It 

II 

IL 



Iwri. . 
(a---r.: 
iaao ( 
:hw n 



u 



or 

^ p 

5 r 
4 c 
•ir 
3 .C 

2c 

2r 

ic 

Ir 



3 . 3 v/ ;C 

ap Tr . 



, la ) 



It 

II 

11 

ti 
II 
1 1 
H 



)V 

5 r 

4c 

4 r 



j 



c 



" 2 r 




I 
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57^! 


Pi'nensionlaau rin^p^r': 




‘37 ' 


para*.ie 'g c r 3 ( ^ i-(7 '’/ 21 - r . ) 


6ror 


57b’ 




5c 


:^8C 




3r 


j‘ j2 




4c 






•i-^’ 






3c 


3'i-. 




3r 


330 




2c 


3^2 




2'^ 


5 




Ic 


530 




1.’ 



53b 


'■'-)r 3o -race '>f A. (initially zero) 


300 


(Zr«^r) for c •■mpu 


radiufj from (A) 


602 


■■'i’-r.iporary 3^'. ’raL^ i’ ;r 
'".;flcc3 )V iat-rfac-:. 


vai'i?.3 a : c r.-.- 


604 


Inooi*:- all /a if cri'- r 


in (A): 6 


’ 'Uu 


0 V ra ye for <jxo ri ^ n o j. 


al wdbr - O^iiic 


00- 1 


0 •'orayo [’ir oc'-iara !’• 


a; oabr . ■ nLi ■- 


61'^ 


Jl 4 1 cT* ' 6 • ^ ^ C 0 Ip *>J. j *3 


. f r .a V i 


. 1 1 


i]s^c: r.j .L .'•r* 


'■.vora ylno 


M , 


! :nus>^u 






X 1 1 i i'. 1. a j. 1 y z i:* r 0 - ■ ■• 3 ' r .a .•• 


, par i. . , Grz 




■f oac 4 :r-oap o irli; 


.r. : n or • "rai.’. 


52 


Inztjr'' 'i . 1.* Z-'’" vj 


>' ’ ’’C - 


'■•2C- j't 


'■-'r.iDoo 




\ 


i.n A' 






Iao_r ;cal: fac3 r- 


V- - . l •. i,; valu"' 




' ccn a r.i i-iO ;.i ■ "i.ia i 


•'I n (.uau . ao i-1 




r' i'l’ 0 .lOi'i'.x i r^-ac 4 
>' jv faz''-} 


'.1 .' xr in ..or":orj 


3 '1 . 


( . .0, ' jf iajor'ac ' -o 


\ 

^ J 



X 
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FLOW SHEET FOR SIX G-ROUPS 
(Superscript i denotes cnercy C^oup, fron 1 to 6) 
(Subscript n denotes space point, fron 1 to 20) 





0 



<D 
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FLOV; 3HEFT FOR SIX SROLFS (Continued) 
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FLOW SHEET FOR SIX GROUPS (Continued) 
SUBROUTINES 
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GHAPTER III 



RESULTS 



I • Re:ictor Corrtiosltlona Studied 

a. Core Materials « In order to simplify the 
problem, a choice was made as to the type of fuel element, 
coolant and moderator. The choice was based on a considered 
adaptability to high output mobile reactors. High power re- 
quirements, coup^-ed with a need for light weight and small 
size, dictated the use of hlp’hly enriched uranium and clad- 
ding materials v/hlch v^ould withstand high, temperatures. 

The fuel element was a dispersal element of 

uranium dioxide in a matrix of stainless steel. Some calcu- 

1 

Itatlons have been made on this tyne of fuel element , and good 

high temperature properties are envisioned. The amount of 

uranium dioxide in the steel matrix is kept low in orde"^ to 

minimize radiation damage and to Improve heat transfer 

properties. A similar fuel element, v/ith the steel replaced 

by aluminum, vms operated at the Genera conference on the 

peaceful uses of atomic energy. 

The uranium considered is enriched to 90 % U-235* 

Tho density^ of UOg is 10.95 gm/cc. 

Stainless steel is also used as the cladllng for 

3 

the fuel element. Its choice is dictated by the fact that 
it will withstand sodium up to about l600°F. Its absorption 
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cross section presents little problem with tl.e fuel enrich- 
ment considered here. 

Berylliuia is chosen as moderator because of its 

Good moderatinG properties, licbt weight, and good thermal 

properties. It is protected from the sodium coolant by stain- 

3 

less steel, since beryllium withstands sodium only to about 
1000 aiid higlier temperatures than this are considered. 

Choice of sodium as the coolant is based on its good 
heat removal characteristics and high boiling point. A high 
boiling point eliminates the need to pressurize the core. 

The reflector is considered , to be composed of 
various proportions of Be, Na, and stainless steel. Reflector 
cooling is necessary since some five per cent of the heat 
may be liberated in the reflector. This amount of heat will 
be considerable for higii output reactors. Stainless steel 
is required, as in the core, to protect the beryllium from 
the corrosive action of higii temperature sodium, and for 
structural purposes. 

With a basic reactor chosen, variations in the 
proportions of the materials were studied in order to 
determ.ino their effect on the power distribution. Table III 
of Appendix E suimnarlzes reactor compositions studied. 

The effect of reactor poisons was not considered, 
nl- of the reactors were assumed to be ho ogeneous and in 
the steady state. 
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b. Cross Sections of Reactor Materials « The 
cross sections of the materials considered are listed in 
Tables I and II of Appendix B • They are based on data 
from “Neutron Cross Sections," BNL 325 (August, 2.955 » 

Geneva Supplement ) • 

Reactor Conf ifcu ration « For reasons of mathe- 
matical simplicity the reactor considered is spherical, 
with a reflector of uniform thickness. Because of angular 
symmetry the spherical reactor may be resolved to a one- 
dimensional problem. Cylindrical or rectangular geometry 
would, on the other hand, require two- or three-dimensional 
analysis. While it may be argued that spherical reactors 
will not find extensive use, the results from an analysis of 
spherical reactors can be extrapolated to cover reactors of 
other configurations. 

The magnitude of the multi-dimensional problem 
is such that slm.pllficatlon to a single dlm.enslcn is made 
frequently. Other than the spherical case considered here, 
the only one-dimensional geometry is that of the infinite 
slab. Spherical geometry Is superior, since it gives 
results which may be applied to reactors of other configu- 
rations whose dimensions along various axes are not greatly 
different. 

Fuel Concentration * The proportion of fuel in 
the reactor has an important effect on the radius of the 
critical reactor. It is also related closely to the power 
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level . 

For the mobile power reactor envlaioned here, 
there are certain restrictions on the size \irhich must be 
considered. It cannot be so large that its mobility is 
questionable, and yet it must be large enough so that 
enough coolant may be pumped through to remove the required 
heat , 

From elementary heat transfer considerations 
(Appendix A ) , a core radius of from two to four feet was 
found to be sufficient for heat removal, and since this 
size appeared reasonable for mobile applications, the amount 
of fuel was adjusted In order to give a radius within these 
dimensions . 

Q* Moderator Proportion , The amount of moderator 
considered varied from zero to sixty per cent of the core 
volume. This had the effect of changing the reactor from an 
essentially thermal one to an intermediate one, in which the 
greater part of the fissions take place at energies appre~ 
clably above thermal, 

f. Reflector Composition . It v;as determined that 
as little as two per cent of the reflector volume was re- 
qxiired to be sodium in order to carry off heat generated 
within the reflector. The amount of stainless steel was 
arbitrarily taken as five per cent. 

By assigning the remainder of the reflector volume 
to the beryllium, it was found that the neutrons returned to 
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the core were largely thermal. In order to decrease the 
number of moderated neutrons reflected back to the core, ad- 
ditional sodium was Introduced into the reflector and the 
beryllium decreased, so that the percentage of beryllium 
considered ranged from a low of 55^ to a high (one run only) 
of 95^. 

• Effect of Various Parameters . 

Fuel Concentration . Apart from its effect on 
the critical radius of a reactor, the fuel concentration 
plays an Important part in determining the energy spectrum of 
neutrons in the reactor. Neutrons are born fast and, during 
the slowing down process, may be absorbed by reactor ma- 
terials. Absorption during slowing down Increases the number 
which must be born fast in order that enough remain after 
slowing down to cause further fissions and perpetuate the 
chain reaction. 

An examination of the cross sections of the ma- 
terials considered here shows that the predominant absorber 
of neutrons having energies greater than thermal is the 
uranium. Increasing the fuel thus increases the flux in 
higher energy groups. 

Actually there are three things which may happen to 
a neutron while slowing down: it may be absorbed; it may be 

lost through leakage; or it may undergo scattering collisions 
and be slowed down. The fuel concentration is so small 
that it can have little effect on the slowing down power. 
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The fuel's Important effect la on absorption tendencies and, 
as a consequence, on the leakage, since an Increase in fast 
flux with more fuel Increases leakage. 

It is Important to note that of the neutrons ab- 
sorbed by the uranium during the slowing dovrn process, 
rouerhly four-fifths cause fission of the absorbing nucleus. 
Thus more fuel has the added effect of changing an originally 
thermal reactor toward a reactor in v^hlch fissions at 
energies above thermal predominate. 

Figure 1 demonstrates the above effect of fuel con- 
centration on the ratio of fast to thermal flux. The ratio, 
calculated at the reactor center, is plotted versus fuel 
concentration for twenty different reactors all having at 
least 30^ beryllium in the core. 

It should be pointed out that the ratio of fast to 
slow flux has an Important effect on the ability of the 
reflector to assist in levelling the power distribution. A 
greater number of fast neutrons leaking into the reflector 
Increases the number of slowed neutrons which the reflector 
can return to the core to compensate for slow leakage , which 
would otherwise depress the pov^er at the edge of the core. 

Moderator Concentration . As mentioned above, 
the ratio of the fast to the slow flux depends on leakage, 
absorption, and slox^rlng down. Changing the moderator pro- 
portion in the core affects the first and last terms. For 
simplicity of discussion its effect on leakage may be 
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nei^lected, If the central regions of the reactor are con- 
sidered. 

In comparing the remaining two terms — absorption 
and slowing dov/n--it is importint to note that the criterion 
for cou.parison is the ratio of either to the sum of the tv/o. 
That is, the total probability of one or the other hapo>^ning 
(excluding leakage) must be unity. For this reason, de- 
creasing one has the effect of increasing the other. 

The beryllium moderator is the predominant con- 
tributor to t slowing dov/n term by virtue of its high 
value of ^ . Any decrease in beryllium thus tends to en- 
hance the absorption lii uranium during slowing down. Thus 
the ratio of fast to slow flux varies inversely with the 
berylll^u?. concentration as shovm in Figure 2 where 5^ast/ ^h 
is plotted as a function of per cent beryllium for reactors 
of various compositions. 

c . Zffect of fuel and moderator concentrations 
On lov/er Distribution 

As discussed above, either more fuel or less 
beryllium increases the ratio of fast to slow flux and thus 
raises the ability of thf reflector to, return slov/ed neutrons 
to the core and assist in levelling the power distribution. 
Fig. 3 shov/s this effect for various boryllivim concentrations 
in an otherwise identical reactor. aS the beryllium decreases, 
the pov/er near the core-reflector boundary rises. 
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The uniforT.lfcy ox’ pov'er dletribixtior., measured by 

the rtitlo of maxlriun to avei*:ge power, depexidx. or the point 

at v'hich P occtirs» This may be either at tne center, rhe 

edge, or, iaeally, oo+h simAu taneously , That is, i should 

be adjusted to e pi.ri the center pow^r 3r or'dei'’ to increase 

P '■s Piuc-h ef5 pooeiMe without rf rnui t'-'.neouGlv raising P 
^.vg ' ^ in^x 

iis an ixxai cation of this, the ratio of P /P in Pi>’. 3 

' max' avg ^ 

decreases to l.l6(as berylliun df^creases from JO to 20 per 
cerit), t>xen rises to 1.62 as ’beryllium is fnrthei- decreased. 
It\ tVie latter case th<=^ edge oowor is limltifu,. 

Fig. ^ snows a similar effect on thr po>?er dis- 
tribution for various fuel concentrations. Figs. 4a and 4b 
indicate flux variations for vax'ious fuel concontratioios , 

Fig. 4b s’no" 'S tcat with r ore fuel the greater fast flux 
results in f, significant number of neutronr bein'" returned 
to the core before beiiu co'''pletoiy moderated. These aid 



in levellirs, the power J.istri bution. 

Also it- rig, (^‘.5n fuel) the l^rge ratio of 
fast to thermsl flux results in some 16,t of fissions 
taking piece in the two hic;hor energy i^:rnups, as coxitrasted 
v/lth OTuly 9/-’ the reactor of Fig. 4a (u.d'o fueli, 

d . Hef I ector Comros i ti on 

The s-ffact of v'''rylng +"he proportions of the three 
reflector oomponrn*'« x-r?r invert igp ted. In genor* ’ , for* the re 
fls'ctor composi tions studied , bsryillum makas the l^'rgor.t 
contribution to the slowing io n terra and ahe smallest to 
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the absorption. Therefore, decreasing the reflector's 
beryllium concentration effectively Increases the absorp- 
tion tendencies while decreasing the moderation. Reflectors 
v;lth less beryllium thus return fewer slowed neutrons to 
the core . 

Fig. 5, which shows the power distributions for 
three reactors with identical cores but with differeht 
proportions of berylliu:’! in the reflector, demonstrates 
this effect. The rise in power occasioned by neutrons 
returned from the reflector is greatest for the reactor 
having the most beryllium, in the ref]ector, least for the 
reactor with least reflector beryllium., 

e. nef lector Thickness 

Figs. 6, 7 , and 8 show the spatial distribution 
of flux for reactors with reflector thicknesses of 10, 20, 
and 50 cm. respectively. Since fast neutrons v/hich leak into 
the reflector are subject to a greater chance of slov/ing 
down in a thick reflector, the therm.al flux rises more 
in such a reflector. In turn, this causes a power rise near 
the core-reflector interface, since more therm.al neutrons 
are returned to the core. 

This provides a means of increasing the thermal 
flux in the outer core region where the pov/er normally is 
less than at the reactor center. By a proper clioice of 
reflector this increased thermial flux in the outer core 
region can result in a pov/er level as high as the power 
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at the reactor center. This represents the closest possible 
approach to a unifonr. power distribution for a given homo- 
geneous reactor, since the power distribution sags between 
these tv/o maximums. 

Fig. 9 shows the effect of varying the reflector 
thiclcness on the power distribution for a given core, A 10 
cm reflector gives a power peak at the reactor center. 
Somewhere between 10 and 20 cm reflector thickness the 
optimum pov/er distribution for this reactor occurs. For 20 
and 30 cm reflector thicknesses a power peak occurs at 
the surface of the spherical core. V/ith Increasing reflec- 
tor thickness a limiting value of Pjnax/^avg is 
approached asymptotically. A plot of the neutron flux in 
Fig, 8 for the ^ cm reflector shows that the flux in all 
energy groups approaches zero asymptotically to the axis. 
Since the curx’ent density J = -Dd(NV)/dr (as given in 
Chapter I), a nearly zero slope to the flux curve means 
that there are very few neutrons escaping from the outer 
reflector surface. This Indicates that a further increase 
in reflector thickness would have little effect on the 
pov/er distribution, so that this is essentially an 
infinite reflector. 

The positive slope (at the core-reflector inter- 
face) to the thermal flux curve means that J is negative, 
indicating a net flow of neutrons from the reflector into 
the core, it is noted that with reflectors of 10 and 20 
cm thickness there is a net leakage of group III neutrons 
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into the reflector. However, with a 50 era thick reflector 
there is a positive slope to the group III flux at the 
Interface, indicating a net return of group III neutrons 
from the reflector to the core. This tends further to in- 
crease the power rise in the outer core regions. 

The significance of fissions in group III can be 
farther Investigated. Figs. Aa,4b, 6, 7, and 8 show flux 
plots for different reactors. Histograms in these figures 
shov; that approximately 60^ of the total fissions take 
place in energy group III. Thus, one would expect that the 
spatial power distribution would follow the core flux plot 
of group III quite closely. An examination of the figures 
mentioned shows that the power distribution does follow 
the group III flux very closely except in the outer core 
regions where the thermal flux begins to rise. At this 
point the power curve also exhibits an upward trend. This 
reflects the fact that the thermal neutron fission cross 
section is approximately three and a half times that of 
the average group III cross section, so that a small inc- 
rease in thenaal flux has great significance. 

Another trend which can be noted from the figures 
cited is that the group III flux distribution sags con- 
tinually as distance from the center increases. Thus 
po’wer which depends primarily on group III will sag also, 
even though the thermal flux remains relatively flat. If 
more moderator is added so that the reactor approaches 
more closely a thermal reactor, the sag in the thermal 
flux in the outer regions will not be compensated for 
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by a return of reflector moderated neutrons, since, with 
relatively few fast neutrons present there will be little 
fast leakage input to the reflector. Fig. 4a Illustrates 
this trend. 



f . Effect of Varying Na/Fe Ratio 
In Fig. 11 is plotted power distributions for 
three reactors identical except for the sodium and stainless 
steel content of the core. The curves are very similar 
except for the outer core region. 

Because Fe has a higher absorption cross section 
it is to be expected that increasing the core Fe content 
will tend to raise the fast flux sornov/hat, thereby 
Increasing fast leakage and the number of slowed neutrons 
returned from the reflector. This point is supported by 
the plot of run (11) v/hich has the highest value of power 
at the corc-r fleeter interface, and the highest Fe con- 
tent. 

An additional effect is the variation in diffusion 
coefficient with varying Na/Fe ratio in the core. Since Fe 
has a higher scattering cross section (by about a factor 
of three) than sodium, and a hlglier atomic density (by a 
sixmllar factor) a greater Fe content tends to decrease D 
and thereby decrease leakage. It is thought that this 
partially explains the sllglit difference betvjeen the plots 
of runs (1) and (10) as shown in Fig. 11. Run (1), with 
less iron and more leakage experiences a slightly greater 
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rlse near the corc-ref lector boundary. 

Since the two effects mentioned are in opposi- 
tion the above discussion is rather hypothetical. In 
general, varying the Ila/Fe ratio has little effect on the 
overall power distribution. 

6 • Optimum Power Distribution 

In Fig. 10 are plotted two optinum pov/er distri- 
butions obtained. Shapes of the curves, even though quite 
different core compositions were considere , are practi- 
cally Identical. .This tends to indicate that a limiting 
value of pow'er uniformity has been met. In other v/ords, the 
power sag observed evidentally cannot be further lessened 
for the reactor type and size considered. 

The tv;o reactors considered in Fig. 10 have 
critical radii of 42 cm and 60 cm. Since reactor size 
affects the proportion of the core which may be affected 
by reflector moderated neutrons, it is thought that size 
will play an important part in determining uniformity of 
power. A reactor with a smaller core v/ill have a larger 
proportion of vo^uine which experiences a power rise due to 
reflector moderated neutrons. However it can bo seen in 
Fig. 10 that decreasing the size from 60 to 42 cm has 
essentially no effect on power distribution. 

The two reactors considered i.i Fig. 10 have 
different material compositions and different critical 
radii. They are representative of the class of reactors 
investigated. It can be concluded, therefore, that in 
the range of size and composition investigated the optinum 
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power curve has a characteristic shape which is that 
ghown in Fig. 10. The ratio of Pmax/^avg this curve is 
approximately 1.16. 

h . Comparison of Four and Six Q-roup Methods 

Fig. 13 gives a comparison of calculations by 
four group and six group methods for the same reactor. 

In genex-’al the results are in good agreement, since the 
critical radii differ by loss than four per cent. The ratio 
of fast to thermal flux was for four groups and 4.7 
for six groups. These values are of the same order of magni' 
tude. 

The plot of spatial pov/er distribution in the 
core shown in Fig. 13 gives different values of Pmax/^avg 
for the two methods, although in general the curves have 
similar characteristics. Better agreement in pov;er dis- 
tribution between methods seemed to exist for reactors 
having other compositions, although other runs on identi- 
cal reactors were not made due to time limitations. Since 
both methods use tlxe age diffusion equation (with minor 
differences v/hich are covered in the discussion of results) 
the results should be very similar. Experimental results 
for I'la, 3 e, Fe -intermediate reactors v;ere not available 
for a Chech on the accuracy of the results. 
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DISCUSSION OF RESULTS 



The principal results of this study have been tiie 
detcriniiiation of critical radius, neutron flux and fission 
pov;er distribution, and the effect upon these quantities of 
varyinG the material compositions of the core and reflector 
as well as varylns the reflector thickness. 

Tue meti.od of attack involved r.odifyinG the a^e 
diffusion equation so as to permit solutions to be ob- 
tained by use of Whirlwind I dicital computer. One method 
v/as developed in Chapter I which used a finite difference 
approximation to tiie diffusion equation v/ith neutrons 
divided into four enerGy Groups. A second method was devel- 
oped in Chapter II which used a method of numerical inte- 
Sration with neutrons divided into six energy 

the results of calculations based upon identical 
reactors using the two different methods were in substantial 
asreement. however, minor differences in spatial power dis- 
tributions and critical radius were noted for solutions by 
the two methods. 

The explanation for differences in the radius and 
power distribution given by the two methods is thought to 
be the result of differences in calculated average cross 
sections for the various energy groups. The energy groups 
are so wide that it is necessary to obtain an average cross 
s<=ction and diffusion coefficient for each energy group. The 
process of averaging could result in incompatibilities 



between average material properties calculated for the 
two methods, particularly in regions v/here the cross sec- 
tion changes rapidly with energy, as is the case for the 
scattering resonance of sodiam near 2700 ev. 

Another major difference between the methods 
involved the determination of the elastic scattering 
input from energy group i to the next lov/er energy 
group i-*-l. In the four-group method a factor f was chosen 
such that f The solution was not sensitive to 

tne value chosen for f as shown by successive runs using 
improved values of f. Chapter I (page 4) discusses the 
manner in which f can be improved from run to run on a 
trlax and error basis. The use of the f factor gives good 
results for v/lde ener y groups. The disadvantage of 
having to estimate f can be overcome by considering a 
greater number of energy groups where the variation of 
the slowing dovm density, q, over the group can be consi- 
dered to be linear. In this case sufficient accuracy 

1 , in out. 

could be obtained by the ass’omptlon that ~ tCQi ) • 

avg ^ 

In the six group method, hov/ever, the assumption 

is that the elastic scattering out of an energy group is 

equal to the average value of the slov/ing dovm density 

in an energy group, This is eouivalent to 

■^avg 1 

assuming that the f factor ecuals unity. These two 
different assumptions concerning the elastic scattering 
out of an energy group probably account principally for 
any differences betv/een results obtained by the two 



methods. This difference can be illustrated by the fol- 
lowing figure: 




The solid line represents the true distribution of slov;ing 
down density. The slowing down density, q, increases for low 
lethargy values corresponding to the birth energies (10 to 
0.1 Ilev) of fission neutrons. For lethargies greater than 
the birth lethargy, q decreases due to absorption and 
leakage. For the relatively wide lethargy groups shown it 
can be seen that the assumption used with the six group method 
that qf^® = is not nearly as true as the representation 

of the true continuous slowing down curve by the assumption 
that qf^^ = f q^*^^ with a value of f chosen so as to 
approximate the true curve as closely as possible. 

The disadvantages of using a relatively small 
number of energy groups are as follows: (1) the true 
elastic scattering out of an energy group is difficult to 
estimate coriectly; (2) determining average neutron cross 
sections for a particular lethargy group introduces errors, 
particularly when cross sections vary rapidly with energy. 
Considering a larger nu.ber of groups greatly reduces 
these troblems. 

The following advantages result from using a 



-S7- 



sinal-L nuir.ber of energy groups: (1) it is possible to 
neglect elastic scattering across groups (for the f'-aterials 
considered); (2) the energy spectrum of fission neutrons 
need not be known v;ith great accuracy if group widths are 
such that all fission neutrons can be considered to be born 
in one or two energy groups. 

.he four group program is arranged to include the 
effects of inelastic scattering. Kov/ever, in this inves- 
tigation inelastic scattering is neglected, since it is 
considered to be of only minor importance for group widths 
and materials considered. 

The plots of neutron flux and power distribution 
gave consistent results indicating qualitative trends as 
parameters v/ere varied. The parameters considered v/ere; 

(1) fuel concentration; (2) reflector composition; (3) Ka/Be 
ratio; (4) ila/pr ratio; and (5) reflector thickness. An 
effort Was made to choose reactors which would have the 
maximum power occurring equally at the reactor edge and 
center to get the smallest value of Pmax/^’avg (ideally 
equal to unity) . The best value that could be obtained was 
1.16 for run ( o) . Fig. 11 shov;s plots for two optimum 
power distributions. Comparison of these tv/o curves with 
other power distributions near the optimum seems to indi- 
cate that the general shape of the curve cannot be cha.nged 

very much for the range of parameters investigated. Hov;ever, 
these results would not apply to reactors having a much larger or 
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smaller critical radius since the shape of the optlraun 

power cui-ve would very likely he considerably chancod with marked 

chanses in reactor size. 

, A ratio of 1.16 for r’max/^avg best value 

obtained for the reactor size and type discussed. 

The two programs which were developed for use v/ith 
the Whirlwind I digital computer operated very satisfactorily, 
and it is hoped that other Investigators will find these 
workable programs to be helpful. Both programs were character- 
ized by rapid convergence to the ansv/er, averaging only 
about four or five Iterations before the convergence criteria 
were met. 



For further work on the subject of uniform pov;er 
distribution, it is suggested tnat multi-region reactors be 
investigated. The Vfhirlwind programs developed In Chapters I 
and II can be modified so as to apply to reactors liaving 
several regions of uniform material properties. The programs 
would lengthen considerably, however, and it probably would 
be necessary to use auxiliary Whirlwind storage. 

It would be of interest to compare the results of 
this theoretical investigation with experimental results in 
order to determine how closely a finite difference or 
numerical integration approximation to the age diffusion 
equation will represent reality. 

An extension of this study to investigate reactors 



of a different range In size using other moderators 
and coolants would supplement the results given here. 
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SUKL^JlRY 

It was the objective of this thesis to study 
the practicability of achlevins a unifora power distri- 
bution in a homogeneous nuclear reactor, through return 
of moderated neutrons from a reflector to the core to 
compensate for leakage, I-iultl-group methods of solution 
of neutron diffusion equations were developed. While it 
was impossible to compare the theoretical results obtained 
with experimantal data, because of lack of such data, the 
results appear reasonable in that they demonstrate trends 
which are in agreement with physical considerations. 

Studies were based on a single basic reactor 
chosen because of suitability to high output mobile 
applications. 

It was found, for the reactor type and size studied 
that, wnile complete unlfonrnity of power could not be 
achieved, nevertheless it was possible to approach 
uniformity quite closely. The best value obtained for 

^raax/^avg 1.16, 
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APPSNDIX A 

For a typical power reactor, Glas?tone^ gives an 
estimated volumetric power output for the core region of 

3 

600 watts/cm. • will be of interest to determine the ap- 
proximate power potential of the rea.ctor of Run 8, ’-^hlch 
gave optimum power distribution with the maximum power 
points equal at the core edge and center. The critical core 
radius was 6O.OO cm. with sodium comprising per cent of 
the core volume. 

Because of time limitations, detailed heat transfer 
calculations are not made. However, an estimate of the 
power potential of the reactor can be made based upon the 
follov/lng calculations j 

Core radius - 60.00 cm. ^ 1.97 ft. 

Core volume - 9 *05 x 10^ cm^ 

Total power output at an assumed 600 watts/cm^ 

= 543 MV 

It will be assumed that the sodliam coolant speed is 

5 ft. /sec. For a mean temperature of 800°F , the sodium has 

2 

the follovring properties : 

^ Z 52 Ib/ft^ 
c =0.31 Btu/lb°F 

1 — S, Qlasstone, “Principles of Nuclear Reactor Engineering," 
D. Van Nostrand Co., Inc., Nev/ York, 1956 , page 636. 

2 — "Reactor Handbook: Engineering," U. S, Atomic Energy 

Commission, August, 1955 > (Geneva), page 266. 
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For pov/er applications, it is hypothesized that a 
cylindrical reactor of 1.97 ft. radius can be chosen which 
will have a similar power distribution and volume to the 
spherical reactor. For 5^3 power output, the mean tem- 
perature rise, At, of the coolant flov;lng axially through 
70 ^^ of the reactor cross section is: 



T = 



L Btu 

5^-3 MW X 3.^13 X lO*" hFI^ 



5 55 ^ == 3600 



Y t^’> 



lb. 






= 546°F 

Thus a power outpiit of 5^3 MW can be obtained from 
this reactor if sufficient heat transfer surface exists to 
give a mean temperature rise of 5^6°F with a sodium coolant 
velocity of 5 ft/sec. With an inlet temperature of 600^F 
and a mean outlet temperature of 1150 *^F , the maximum tem- 
perature of a thin fuel element should not be excessive* 

The ability to transfer tho required heat from the 
fuel elements to tho coolant is the primary limiting factor 
which must be studied in determining the maximum pov/er output. 

It is of Interest to compare the output of this re- 
actor with the power requirements of a B-52 t^/pe Jet boiaber. 

At 563 miles per hour, eight Jets of 10,000 lbs. thrust each 
develop a total of 120,000 thrust horsepower. If an overall 
efficiency of 20^ is assumed, this is equivalent to a heat 
source developing 450 MV/« This can be compared to the 
roughly estimated reactor power output of 5^3 MW. 
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TABLE I 

MIGROGCOPIC CROSS SFGTIORS^" 



Scattering Cross Sections, cTs 

Energy Lethargy, u= 

E (Mev) In (10/E) Ka 


(bams) 

Fe 


Be 


U-235 


0 


10 




0 


2.4 


3.0 


1.7 


6.55 


1.25 


5.62 




0 . 576 


2.4 


3.6 


1.0 


6.45 


1.5 


3.16 




1.151 


2.4 


3.3 


2.1 


6.45 


1.5 


1.776 




1.727 


2.5 


2.9 


2.5 


5.24 


1.7 


1 




2.303 


3.5 


2.8 


2.9 


5.08 


4.0 


0.562 




2.^79 


4.0 


3.0 


3.5 


5.78 


3.5 


0.316 




3.454 


3.2 


3.0 


4.0 


7.24 


4.0 


0.177O 




4.030 


4.0 


3.5 


4.9 


8.7 


3.5 


0.1 




^4.606 


3.2 


4.0 


5.7 


9.76 


3.5 


60.01 X 


10" 


■^5.123 


4.7 


3.6 




10.2 


3.6 


36.21 


\\ 


5.640 


4.1 


6.4 


5.9 


10.7 


3.7 


20.90 


tl 


6.157 


4.2 


2.4 




10.9 


3.0 


12.70 


U 


6.674 


4.9 


4.7 


5.9 


10.9 


3.8 


7.59 


It 


7.191 


7.7 


8.8 




11.1 


3.8 


4.55 


tt 


7. 708 


17.5 


6.2 


6.0 


10.7 


3.8 


2.69 


II 


C.225 


200.0 


6.8 




10.9 


3.8 


1.65 


II 


8.742 


10.0 


7.8 


6.0 


11.7 


3.8 


0.955 


II 


9.259 


^.5 


9.3 




14.6 


3.8 


0 . 569 


II 


>.776 


3.2 


10.5 


6.0 


18.2 


3.3 


0.346 


II 


10 . 293 


3.1 


11.0 




20.6 


3.8 


0 . 2oG 


II 


10 . 610 


3.1 


11.4 


6.0 


19.2 


3.8 


55.0 X 10"*^ 


12.142 


3.1 


11.4 


6.0 


17.6 


3.8 


14.45 


II 


13.474 


3.1 


1-.4 


6.0 


19.0 


3.3 


4.07 


II 


14.806 


3.2 


11.4 


6.0 


19.0 


3.8 


1.00 


II 


16.138 


3.2 


11 . 4 


6.0 


3.6 


3.8 


0.269 


n 


17 . 470 


3.2 


11.4 


5.9 


2.5 


3.8 


0.070 


It 


18.800 


3.2 


11.4 


5.9 


10.0 


3.8 



* lieutron Cross Sections . BIIL -3255 U.S. Atomic Energy 
Com...lsslon, Aug. 1955 (Geneva Supplement) 
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::iCROSGO:iG cross SECTIOrS (continued) 



3. Absorption and Mission Cross Sections, cr ' g ^ and cTf (barns) 



Energy Lethargy, u* U-235 



E ^ _(i’:ev) 


ln(10/E) 


I'a 


Fe 




OTf 


cr & 


0 


10 


0 


0 


0 


0 


1.29 


1.55 


0 


5.62 


0.576 








1.29 


1.55 




3.16 


1.151 


0 


0 


0 


1.29 


1.55 


0 


1.778 


1.727 








1.30 


1.56 




1.000 


2.303 


0 


0 


0 


1.27 


1.52 


0 


0.562 


2.879 








1.20 


1.44 




0.316 


3.454 


0 


0 


0 


1.30 


1.56 


0 


0.177S 


4.030 








1.50 


1.80 




0.1 


4.606 


0 


0 


0 


1.70 


2.04 


0 


60.01 X 10- 


5.123 








2.00 


2.40 




36.21 ;; 


5.640 


0 


0 


0 


2.35 


2.82 


0 


20.90 " 


6.157 








2.80 


3.36 




12.70 " 


6.674 


0 


0 


0 


3.4 


4.08 


0 


7.59 


7.191 








4.1 


4.92 




4.55 


7.708 


0 


0 


0 


5.0 


6.0 


0 


2.69 


8.225 








6.0 


7.2 




1.65 - 


8.742 


0 


0 


0 


7.4 


3.9 


0 


0.955 " 


9.259 








9.1 


10.9 




0.569 


9.776 


0 


0 


0 


11.5 


13.8 


0 


0.346 " 


10.293 








14.5 


17.4 




0.200 " 


10.310 


0 


0 


0 


19.0 


22.8 


0 


55.0 X 10"^ 


‘ 12.142 








52.0 


62.4 




14.45 " 


13 . ^-^74 


0 


0 


0 


55.0 


66 .0 


0 


4.07 " 


14.806 








30.0 


36.0 




1.00 ” 


16.133 


0.1 


0.1 


c 


72.0 


86.4 


0 


0.269 


17 . 470 


0.15 


0.1 


0 


150.0 


187.5 




.0.070 ” 


18.800 


0.30 


0.6 


0.006 


320 : 


375 


0.0001 
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TAPLE II 

AVEkaGE hICRCSCOFIC CROSS SECTIONS 
PER EilERGY GROUP 



■» 



A. Four Energy Groups (cross sections in bams) 







^ tr 


c^a 


cFf 


J/ 


5 


Group I 


( Lethargy 


range 0 


to 4.606) 








Na 


3.10 


3.01 


0 


0 




0.0845 


3e 


3.175 


2.94 


0 


0 




0.209 


Fe 


3.20 


3.16 


0 


0 




0.0359 


0 


2.76 


2.64 


0 


0 




0.120 


b-5 


6.6^1- 


6.64 


1.597 


1.331 


2.46 


O.OO 85 


Group II 


( Lethargy 


range 4 


.606 to 10.81) 






i^ia 


22.22 


21.6 


0 


0 




0.0845 


Be 


5.94 


5.50 


0 


0 




0.209 


Fe 


7.10 


7.01 


0 


0 




0.0359 


0 


3.76 


3.60 


0 


0 




0.120 


u -5 


12.90 


12.90 


7.65 


6.54 


2.46 


0.0085 


Group III (Lethargy range 


10.81 to 


13.80) 






::a 


3.16 


3.07 


0.067 


0 




0.0345 


Be 


5.97 


5.54 


0 


0 




0.209 


Fe 


11.4 


11.26 


0.083 


0 




0.0359 


0 


3.8 


3.64 


0 


0 




0.120 


U-5 


13.0 


13.0 


106.2 


£8.1 


2.46 




Group T 


( Lethargy 


18. SO, E 


= 0.070 


ev, 1000 


1 F thermal energy) 


I!a 


5.2 


3.11 


0.30 


0 




0.0845 


Be 


5.9 


5.46 


C.006 


. 0 




0.209 


Fe 


11.4 


11.26 


C.60 


0 




0.0359 


U -5 


10.0 


10.0 


575. 


320 . 


2.46 


0.0065 



Averaged with respect to lethargy. 
U-5 stands for U 255 
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CnOSb SECTIONS (Continuedj 
B. 3ix Groups 



6 ( fast ) ; 


0.234 


r;rv 


to 


10 


j^lev 


5 


5.5 


kev 


to 


0.234 


i^ev 


4 : 


128 


ev 


to 


5.5 


kev 


3 : 


3 


ev 


to 


128 


ev 


2 ; 


.07 


ev 


to 


3 


ev 


Kth) : 


.07 


ev 









Na atomic densit ’’' ( xlo'^^j 

Be 

i^e 

0 (in UO 2 ) 

U-i 35 (QQ-^jin UOo) 

U-23B 



.0254 Ao 


.029 ^ 


.12 


.074 


.0848 


.0119 


.01,01, 


. 0/.16 


.02224 


. 0 J283 


. as 251 


.0028 



.ut '45 

.2J9 

.0359 
.12 
.JJS 5 
. O004 



Absorption Cross-jections ( barns )( includes lission absorption ) 



finer, <yy Group 


Na 


Be 


Fe 


0 


U -235 


U -238 


1 


0.3 


0.0 


0.6 


0. 0 


374.0 


1.6 


2 


0.1 


vJ • J 


0.2 


0.0 


155*0 


0.8 


3 


0.0 


0 • 0 


0. 0 


0.0 


54.8 


0.5 


u 


0*0 


0.0 




0. 0 


12 . 4 


0.5 


5 


0* J 


0. J 


0.0 


0. 0 


3.1 


0.5 


6 


0. J 


0. J 


0.0 


0.0 


1.6 


0.5 


Elastic Ocatte 


rin;^ orooo- 


actions (barns): 




1 


3.3 


5 .C 


11.4 


3.^ 


11.0 


8.4 


2 


3.3 


5 .^' 


11.4 


3.8 


10.5 


^.8 


3 


3.1 


( .0 


11.4 


3.8 


8.5 


^.0 


4 


27.4 


0. 0 


6.4 


3.8 


8.3 


‘ .0 


5 


3.6 


5.1 


2.3 


3.3 


8 .0 


^.0 


6 


2.7 


2 . A 


2.7 


2 . b 


^ . J 


‘^.0 



rission Crosi-iections (barns); 



1 

2 

3 

L 

5 

6 



3-iO. 'J 
12 C .2 

45.7 

IU.4 

2.6 

1.3 
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CROSS SECTIONS (Continued) 

Since BML-325 lists total cross sections for 
elements other than fissionable materials certain 
assumptions are necessary in order to compute scattering 
and absorption cross sections. The following assumptions 
were made with respect to all elements other than 
uranium’; 

(1) the thermal absorption cross sections listed 
in BKL-325 for 2200 a/s neutrons was decreased 
by a l/v relation to the .07 ev thermal 
energy considered here; 

(2) the scattering cross section was assumed to 

be the total cross section at the point 
v/here the absorption cross section, scaled 
by l/v became negligible. , 

These assumptions resulted in absorption in only the 
thermal and next higher groups. 

For the uranium the fission cross section is listed 
in BNL-325* The absorption cross section was synthesized 
from this under the assumption that the capture to fission 
ratio is constant. The scattoi’lng cross section of U-235, 
and all cross sections of U-238, proved to be unimportant 
for this probler . The capture to fission ratio of U-235 
v/as taken to be 0.164. 
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TABIii: HE 

REACTOR COlTOSI'xIOrS 



Volunie Per Cent 

Core Reflector 



Run Ro. 


Ka 


Be 


re 


UOp 


Wa 


Be 


Fe 


1 


60 


30 


9.7 


0.3 


20 


75 


5 


2 


60 


30 


9.6 


0.4 


20 


75 


cr, 


3 


60 


30 


9.5 


0.5 


20 


75 


5 


4 


70 


20 


9.6 


0.4 


10 


85 


c 




70 


20 


9.6 


0.4 


20 


75 


5 


6 


70 


20 


9.6 


0.4 


40 


55 


5 


7 


60 


30 


9.7 


0.3 


10 


35 


5 


8 


70 


20 


9.7 


0.3 


10 


85 


5 


9 


80 


10 


9.7 


0.3 


10 


85 


5 


10 


50 


30 


19.7 


0.3 


20 


75 


5 


11 


40 


30 


29.7 


0.3 


20 


75 




12 


75 


15 


9.8 


0.2 


5 


90 


5 


13 


70 


20 


9.6 


0.4 


15 


80 


5 


14 


60 


30 


9.8 


0.2 


20 


75 


5 


15 


60 


30 


9 . 8 


0.2 


40 


55 


5 


■ 16 


30 


60 


9.8 


0.2 


20 


75 


5 


17 


30 


60 


9.3 


0.2 


40 


55 


5 


18 


90 


0 


9.7 


0.3 


25 


70 


5 


19 


60 


30 


9.7 


0.3 


2 ^ 


95 


2.5 
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TABLE 32 : 



WHIRLV.'IMD INrUT D^Ta 
FOR FOUR ENERGY GROUPS 



Run II 


Enersy 
o . Group 


D 


f 


CQ 






Kc)** 


I 


1.867 


0.756 


0.02825 


0.000106 


0.000219 




II 


0 . 5686 


1.412 


0.07548 


0.000523 


0.00107 




III 


0.9769 


1.44 


0.05243 


0.008781 


0.01444 




T 


0.9868 


1.00 


0 . 5196 


0.03471 


0.05247 


1 (R) 


I 


1.138 


0.68 


0.06154 


0.0 






II 


0.531 


1.08 


0.12218 


0.0 






III 


0.593 


1.05 


0.11539 


0.000692 






T 


0.601 


1.00 


0.10581 


0.00566 




2 (C) 


I 


1.871 


0.50 


0.02819 


0.0001412 


0.000291 




II 


0.569 


1.75 


0.07540 


0.000697 


0.00143 




III 


0.983 


1.28 


0.05233 


0.011136 


0.01929 




T 


0.990 


1.00 


0.05186 


0.04298 


0.06939 


2 (R) 


I 


1.138 


0.56 


O.O6I54 


0.0 






II 


0.526 


1.36 


0 . 12236 


0.0 






III 


0.593 


1.18 


0.11539 


0.000692 






T 


0.601 


1.00 


O.IO58I 


0.00460 




3 (C) 


I 


1.870 


0.50 


0.02627 


0.0001768 0.000364 




II 


0.567 


2.50 


0.07548 


0 .0^0871 


0.001790 




III 


0.982 


1.57 


0.05242 


0.013490 


0.02412 




T 


0.988 


1.00 


0.05194 


0.051228 


0.08760 


3 (R) 


I 


1.138 


0.56 


0.06154 


0.0 






II 


0 . 531 


1.15 


0.12218 


0.0 






III 


0.593 


1.05 


0.11539 


0.000692 






T 


0.601 


1.00 


0.10581 


0.00566 




4 ( 0 ) 


I 


2.215 


0.621 


0.21473 


0.0001412 


0.000292 




II 


0.580 


1.526 


0.065327 


0.000697 


0.00143 




III 


1.189 


1.660 


0.03 810 5 


0.011306 


0.01929 




T 


1.194 


1.00 


0 . 037^16 


0.043668 


0.06989 


■^(c) 


stands for reactor 


core 


properties. 


(R) for reflector 



properties . 
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WHIRLWIND INPUT DATA 
(continued) 



Run 

No. 


Energy 

G-roup 


D 


f 








UR 


I 


1.039 


0.68 


s 

0.06884 


0 






II 


0.522 


1.08 


0.13231 


0 






III 


0.537 


1.05 


0.12968 


0.000522 






T 


0 . 5 ^^ 


1.00 


0.12820 


0.003914 




56 


I 


2.215 


0.44 


0.021473 


0.0001412 


0.000292 




II 


0.580 


2.50 


0.065327 


0.000697 


0.00143 




III 


1.189 


1.57 


0.038105 


0.011306 


0.01929 




T 


1.19^ 


1.00 


0.037816 


0.04367 


0.06989 


5 R 


I 


1.138 


0.57 


0.06154 


0 






II 


0.531 


2.50 


0.12218 


0 






III 


0.593 


2.00 


0.11539 


0.000692 






T 


0.601 


1.00 


0.10581 


0.00566 




60 


I 


2.215 


0.5 


0.021473 


0,0001412 


0.000292 




II 


0.580 


2.0 


0.065327 


0.000697 


0.00143 




III 


1.189 


1.57 


O.O38IO5 


0.011306 


0.01929 




T 


1.19^ 


1.00 


0.037816 


0.043668 


0.06989 


6r 


I 


1.4004 


0.85 


0.04694 


0 






II 


0.5445 


1.22 


0.10209 


0 






III 


0.7501 


1.02 


0.0868 


0.001033 






T 


0.7581 


1.00 


0.08587 


0.006 




70 


I 


1.867 


0.621 


0.02825 


O.OOCIO6 


0.000219 




II 


0.568 


1.528 


0.07548 


0.000523 


0.00107 




III 


0.9797 


1.660 


0.05243 


0.008781 


0,01444 




T 


0.9868 


1.00 


0.05196 


0,03471 


0.05247 


7 R 


I 


1.034 


0 . 68 


0.06944 


0 






II 


0.520 


1.08 


0.1331 


0 






III 


0.535 


1.05 


0.1331 


0.000491 






T 


0.540 


1.00 


0.1293 


0.003773 




80 


I 


2.209 


0.621 


0.02158 


0.000^06 


0.000219 




II 


0.5796 


1.528 


0.06535 


0.0005231 


0.00107 




III 


1.184 


1.660 


0.03813 


0.005871 


0.01444 




m 

J. 


1.190 


1.00 


0.03784 


0.0354 


0.05247 


8 R 


I 


1.039 


0.68 


0.06884 


0 






II 


0.522 


1.08 


0.13231 


0 






III 


0.537 


1.05 


0.12968 


0.000522 






T 


0.544 


1.00 


0.12820 


0.003914 
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whirlwind INPUT DATA 









(conti 


nued) 






Run 


Energy 












No. 


Group 


D 


f 






Jy' ^ f 


9 C 


I 


2.700 


0.621 


0.01436 


0.000106 


0.000219 




II 


0.591 


1.528 


0.05527 


0.000523 


0.00107 




III 


l.i 4-98 


1.660 


0.02386 


0.009121 


0.01444 




T 


1.500 


1.00 


0.02373 


0.03631 


0.05247 


9 R 


I 


1.03^ 


0.68 


0.06944 


0 






II 


0.520 


1.08 


0.1331 


0 






III 


0.535 


1.05 


0.1311 


0.000491 






I 


0.540 


1.00 


0.1293 


0.003773 




IOC 


I 


1.698 


0.9 


0.02919 


0.00011 


0.000219 




II 


0.564 


1.2 


0.07287 


0.000523 


0.001073 




III 


0.779 


1.4 


0.05522 


0.00931 


0.01444 




T 


0.784 


1.0 


0.05474 


0.03904 


0.05247 


lOR 


I 


1.138 


1.0 


0.06154 


0 






II 


0.531 


.95 


0.12218 


0 






III 


0.593 


1.0 


0.11539 


0.000692 






T 


0.601 


1.0 


0.10581 


0.00566 





TABLE V - SUI'IMARY OF RESULTS 
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rH 
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o5 
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1 
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Si 

<D 

O ^ 




MD 


in 
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MO 
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MO 

1 


m 


NO 
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MO 
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AFFEIIDIX C 

Derivation of Equations 

Equation (5»5d.) of Chap. I v,'as stated without proof 
and used as a basis for the boundary conditions of the four 
group method. A proof of equation (3.5d) follov/s: 



Prove that 




As a preliminary step it is necessary to prove the 
following relationship 
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AFPSrPIX C (continued) 



using the relationship 




rt~ f 



Thus, equation (a) is proved. IText introduce the variable ^ 
and use the definition of V in spherical coordinates 






Use equation (b) in (a) with replacing ^ j 
Where f = 



Then 
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Appendlx C ( contl-iued) 
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/I the reoult beco.'nes (continued) 
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hFji'SIIDIX C (continued) 




since is to be identified v;ith ^ ^ . This is the 

proof of equation (3»5h) of Chap. I which is concerned 
only with the first part of the above expression relating 
to evaluation at the point r^^. 
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